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Abstract 


Ka ^h^r.s t i Ca  1 TiOCtls  of  fluid  *rariSrT.issior.  lin^s  zf  arlitrarv  cr'css 
section  were  cevelct-ed  to  describe  the  protagaticn  of  sc.all  signals  ir. 
developed  r.ear.  flow.  The  t;ear.  flow  ynay  oe  either  .^ar.irjar  or  turbulent, 
and  the  fluid  r.ediur.  either'  cor.t  res  sidle  or  incor.pressirle.  Expressions 
for  both  the  frequency  and  tire  dor.ain  cynar.ic  responses  are  derived. 

Tne  analogy  between  electrical  transmission  line  theory  and  transient 
fluid  mechanics  is  used  tc  derive  the  series  impedance,  shunt  ad-r.ittance , 
c.haracter istic  impedance  and  prcpagaticn  cperatcr  for  each  distributed 


parameter  model. 

A characteristic  frequency  and  radius  are  oef 
arbitrary  cross  section  which  relates  the  viscous, 
prcperties  and  the  inertial,  adiaoatic  properties, 
cross  section,  the  characteristic  radius  essential 
circular  line  having  the  same  viscous  to  inertial 
the  ratio  of  the  excitation  frequency  to  the  line 
frequency  and  the  characteristic  radius  used  ir.  dy 
rcughly  ar.alcgous  tc  the  Aeynclds  num.ber  and  hydra 


ined  for  any  line  cf 


isotherm, al  line 


For  an  arbitrary  line 
ly  defines  an  equivalent 
steady  flow  prcperties; 
cnaracteristic 
nam.ic  analysis  are 
ulic  radius  used  in 


steady  flow  ana.ysis. 

In  developing  the  models,  the  frequency  response  of  transm.ission 
lines  cf  annular  cross  section,  not  previously  available,  was  first 
solved  and  comt'ared  tc  existing  solutions  for  circular  and  rectangular 
lines,  as  well  as  experimental  data  s'ubsecuer.tly  generated  for  blocked 
annular  lines.  Secondly,  the  frequency  and  time  domain  solutions 
appiicacle  tc  circular,  rectangular  and  annular  lines  carrying  laminar 
or  no  mean  flow  were  obtained  in  terms  of  equivalent  circular  lines  and 


xvii 


A-:T/DS/r.A.'77:-: 

cor.parec  tc  e>;i£‘.i:.g  thecrv  and  exjer iTT.ental  data.  Finally,  tr.e  ccr.cett 
cf  a characteristic  frequency  and  radius  vas  extended  tc  arbitrary’  lines 
carrying  ccc.rressiale  cr  inccatressicle  turbulent  rtean  flov.  The  fre- 
quency and  tiDe  dor.ain  results  are  shewn  tc  be  coricaratively  ver".’  siir.tle 
tc  artly,  yet  ccrcarable  in  accuracy  tc  existing  vern'  corit-lex  r.odels. 

The  cor.car isens  for  turbulent  flew  were  restricted  to  incer.pressib le 
fluids  in  circular  lines  since  this  was  the  extent  of  the  available 
analytical  and  experir.ental  data  base.  However,  the  mcdel  itself  can  be 
aprlied  tc  noncircular  lines  and  coircressibie  flow.  In  all  cases  the 
r.oaels  are  tjore  versatile  and  easy  tc  apply  than  existing  nodels  without 
sacrificing  accuracy  when  compared  tc  experimental  data.  The  models 


should  therefore  be  of  significant  value  to  the 
fluid  systems  desigrj  and  analysis. 
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-SIGNAL  risponsi  cf  fluil  tkansnissiok  lines 
INCLUDING  DEVELCFED  KEAi:  FLOW  EFFECTS 

I.  Ir-troducticr- 


nc 


The  Grcjagatior.  cf  ST:,all  3igr;als  through  ducts 

lid  has  heen  of  interest  to  scientists  and  engineers  fcr  sorie 
tine.  Witr.  the  advent  cf  fluidics  and  fLuerics  in  the  Late  L?50's  and 
up  tc  tne  present  tine,  nany  investigations  have  been  ccnpleted,  bcth 
anaLyticaily  and  experimentally,  in  order  to  predict  and  measure  the 
parameters  governing  small  signal  propagation. 


T.ne  modeling  cf  fluid-filled  transmission  lines  has  crogressed 
fr-om.  the  basic  lossless  line  model,  tc  the  line  whose  losses  m.ight  be 
called  linear  with  mean  velocity,  to  the  present  freouency-dependent 
friction  models.  Solutions  have  been  obtained  in  the  frequency  do-.ain 
fcr  Lines  cf  circular  arid  rectaiigular  cross  section,  as  these  cross 
sections  are  commonly  found  in  fluid  systerris.  For  the  fluid-filled 
circular  line,  the  transient  response  tc  small  amplitude  impulse  and 
step  inputs  has  been  obtained  in  various  form.s. 

Recently,  investigations  have  been  carried  out  and  models 
developed  to  determine  and  to  describe  the  response  of  liquid- filled 
circular'  lines  to  oscillatory  inputs,  where  the  inputs  are  superimposed 
on  a steady  developed  flow,  and  where  the  steady  developed  flow  may  be 
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either  lat.inar  or  roily  roriulent.  A lir.itei  arr.cur.t  of  exterir.er.tal 
data  .-jaS  troviced  t-ore  decree  of  creder.ce  to  these  tiodels . 

Frcbler  lescrirticr 


Turing  the  literafore  search  rreoeding  this  wc-rV,  it  att-eared 
tnat  a t.a;|ority  of  the  investigations  were  verj"  specialised  ani  that 
an  effort  to  synthesize  the  results  obtained  by  others  was  in  crder. 
huch  of  the  previous  wcrV  relating  to  transient  response  in  the  tire 
dot.ain,  and  tc  ooriined  oscillatcr-y  and  steady  developed  flow,  was 
restricted  tc  liquid-filled  circular  lines. 

A solution  for  the  srall-signal  response  of  fluid  lines  of 
concentric  annular  cross  section  appeared  tc  be  a useful  first  step  in 
synthesizing  and  expanding  the  arrlicabilitp-  cf  transrissicn  line 
modeling.  Since  the  annular  line  includes  the  circular  line  as  a 
special  case  (when  the  ’nside  radius  is  equal  to  zero)  and  approaches 
the  parallel  plate  configuration  in  its  other  ex.treir.e , a solution  was 
first  obtained  for  this  line. 

In  studying  the  solutions  available  for  circular  and  rectangular 
lines  as  well  as  the  solution  for  annular  lines  given  herein,  it  appeared 
that  a Tiethod  cf  expressing  the  frequency  dependent  response  of  lines  cf 
arbitrary  cross  secticn  in  terrr.s  of  only  the  circular  litie  night  be 
developed.  Through  the  use  cf  the  various  analogies  between  electrical 
transnission  line  theory  and  fluid  dynamics,  and  in  particular  by 
constructing  and  analyzing  the  equivalent  electrical  circuit  of  a 
generalized  fluid  transmission  line,  such  a method  was  found. 

It  followed  that  if  the  frequency  dependence  of  a line  of 
arbitrary,-  cross  section  could  be  expressed  in  terms  of  an  equivalent 
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circular  lire,  then:  (1)  Tne  srall  sigr.al  transient  (tire  dcr.ain) 
response  of  lines  cf  arlitrar'/  noncircular  cross  section  right  be  found, 
since  this  response  is  hnown  for  circular  lines;  (2)  The  characteristic 
irtedance  cf  lines  cf  ncncircular  cross  section  ray  be  expressed  in 
terras  of  eq'ui valent  circular  lines;  consequently,  irtedance  ratching  of 
lines  cf  different  cross  section  ray  be  possible;  (3)  The  corputer 
noceling  cf  all  lines  cf  any  cross  section  in  a fluid  syster  right  be 
reducer  tc  roce— mg  tne  ecucva— ent  corcuoai'  ^.nes;  v**)  A retnoc  ct 
deterr.ining  the  srall  signal  response  cf  fluid  lines  of  arbitrary  cross 
section,  carrying  fully  developed  mean  flow  right  be  found. 

All  of  these  possibilities  were  investigated  and  the  associated 
probler.s  were  solved,  subject  to  the  restrictions  specified  herein,  it 
is  nctfed  that  wr.6re  r‘C'S£ihl€  th*-  sirtrl^st  rrjodfel  achi^vcila  is  usst.  ir. 
accord  with  Olcenberger' s observation  that  only  those  models  which  are 
reasonably  uncorpl'icated  find  wide  use.  This  observation  is  particularly 
appropriate  tc  fluid  systems,  which  generally  contain  a myriad  of 
components  arid  lines,  some  of  which  are  usually  ncncircular  in  cross 
section. 


Orgar.ization  atid  Contents 

Since  this  work  covers  a rather  broad  range  of  transmission  line 
problems,  specific  background  information  and  references  will  be 
introduced  as  appropriate.  Further,  where  specific  exa.T.pleE  are  deemed 
necessary',  these  will  be  given  in  terrr.s  of  the  more  common  rectangular 
and  annular  cross  sections  for  which  a relatively  larger  body  of 
technical  reference  material  exists.  Lengthy  mathematical  details  which 
ax'e  of  minor  import  in  understanding  the  various  developments,  as  well 


3 


cc  c^r'tair.  6Xdr.r  l€5  ir.d  cr'ror  sriil'.’soc*  hsv's  *c  th^ 

arr-eridicrE . 

In  Se:tl~r:  II,  the  exact  theor"/  descritir.r  the  sr, all-signal 
^ilu*d  ^T‘3r»STr. * ® ^ * cri  snn’uds.'^  c^i'os^ 

presented.  The  results  are  compared  v.ith  those  oltained  'zy  others  for 
circular  and  rectangular  lines.  Thus,  Section  II  serves  also  as  an 
introduction  to  the  synthesis  of  results  aptlicahle  to  circular  and 
ncncircular  lines.  High  and  lov  frequenc/  apprcxinations  are  developed 
for  the  annular  line,  and  exterinental  results  are  presented  which 
ccnfim.  the  theor-y. 

In  Section  III,  a characteristic  frequency  and  corresponding 
characteristic  radius  are  deterr.ined,  which  allow  the  frequencp*  response 
cf  noncircular  lines  to  be  presented  in  terr.s  cf  equivalent  circular 
lines.  The  results  are  ccr.pared  with  those  previously  obtained  for 
rectangular  and  annular  lines  by  using  the  classical  hydraulic  radius 
together  with  circular  line  theory.  The  results  predicted  using  the 
characteristic  radius  and  circular  line  theory  are  shown  to  be  clearly 
superior  to  the  results  predicted  using  the  classical  hydraulic  radius 
and  circular  line  thecr^.-.  However,  the  hp’draulic  radius  nay  be  ad'fusted 
to  allow  its  use  in  the  circular  line  theerv.  A concert  cf  ir.tedance 


matching  is  also  developed  in  Section  III. 

Many  practical  probleTTis  are  more  concerned  with  the  tiDe  doc.ain 
response  cf  transmission  lines.  In  Section  IV,  the  results  of  others, 
v;hich  describe  the  impulse  and  unit  step  response  of  circular  lines 
only,  are  adapted  to  lines  of  noncircular  cross  section  and  also 
compared  with  existing  experimental  data. 
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Sec-ior.  V is  d^vctec  t:  th*  sir:  lificasior.  ar.d  syrjtr.eEis  cf 
results  dessrilir.g  the  oordinaticr.  of  steady  developed  f^buler.t  rear, 
flow,  ar.d  lar.irar  cscillator^/  flow  produced  iy  a srr.all  siruscidal 
pressure  ir.put  signal.  The  solutions  are  extended  to  include  lines 
of  nonuircular-  cross  section. 

In  Section  VI , the  frequency  domain  results  of  Section  V are 
transposed  into  the  tine  dor.ain,  and  the  various  coriinations  of 
downstreai:.  pressure  and  flow  response  to  upstreat  step  inputs  of 
pressure  :r  flow  are  presented. 

In  Section  VII,  the  last  section  of  this  v:crk,  a sati.T.ar\’  cf  the 
conclusions  and  reconnendations  is  rresented. 


: : . Tr.^  Sr.all-Sigr.al  ?es: cr.se  :f  .».r.r.ular 
."iu:d  Trar.sr dssior.  Lir.es 


Ir.troduCCicr. 

"he  sr.all-sigr,al  frequency  response  of  fluid  trar.sr.iss ior.  line 
of  circular  cross  secticr.  has  bee.',  obtained  in  various  fcrr.s  [1-3]*. 
Kichcls  [3]  expressed  the  results  in  a r.anr.er  sir.ilar  to  tnat  used  for 
electrical  transaission  lines.  Schaedel  [*•]  followed  closely  The 
apprcac.'.  used  by  Kichcls  and  obtained  the  response  for  lines  cf 
rectangular  cross  section.  In  octh  cf  these  studies,  a distributed 
pararceter  line  riocel  was  developed  and  results  were  obtained  for  the 
irjpedance,  adr.ittance , prepagation  operator  and  characteristic  iirpedance 
cf  the  line. 

There  are  a nuir-ber  cf  applications,  such  as  aircraft  pitot-static 
systeriS,  hydraulic  r.ining,  and  instruxentaticn,  where  annular  lines  are 
used  and  the  snail  signal  frequency  response  characteristics  are  cf 
interest.  Tne  purpose  cf  this  section  is  to  present  the  sr.all -signal 
response  of  rigid  pneun.atic  transr.lssion  lines  cf  annular  cross  section, 
figure  1.  The  analytical  approach  taker,  is  sinilar  tc  that  used  by 
Kichcls  [3]  and  Schaedel  The  detailed  analytical  development  is 

presented  in  Appendix  C.  finally,  the  annular  line  results  will  be 
compared  with  these  obtained  for  circular  and  recta;;gulor  lines. 


^Numbers  in  brackets  designate  Bibliography  at  end  of  dissertaticn . 
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Sr.all  Signal  Equationis  and  Snluticn • The  ecuaticns  extressing 
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The  dependent  variailes  are  assumed  to  be  of  the  following  form 

T(r,z,t)  = T + T'  = T + i.T(r,z)  (4) 

u(r,2,t)  = u + u'  = 4 + u(r',z)  (5) 

?(z,t)  =?+:-’  = F + dn(z)  (£) 

wner-e,  for  example,  T is  the  average  temperature  and  T'  is  the  instanta- 
neous temperature  fluctuation.  Sufc'Stitution  of  these  variables  into 
equations  (2)  and  (3)  gives 


a - A'T 


cr" 


I cLZ 
r cr 


-1 


(£) 


Ass'ar.irjg  isctherr.al  walls  and  zero  velocity  at  the  walls,  the  bcundart’ 
conditions  for  the  annular  line  are 


IT  = C,  u = ^ 


AT  = C , u = C at 


(9) 


(11) 


The  DOTT.entur  and  energy  equations  are  thus  identical  in  fort,  anc 
the  solutions  for  the  dirjensionless  velocity  and  ter.perature  profiles 
are  Bessel  functions  of  zero  order. 


:uiu 


1 9 ho 
p 9z 


= 1 + 5 J f 


IF-")  ^ 


(11) 


" ^^^0  1 


u . o 
Vn* 


+ FKj, 


( {5..f 


(12) 


The  arbitrary  coefficients  (2,1,1,  and 
conditions  are 


r)  cotacnec  fro.t.  t.te  ooundar\’ 


Ko 


£ = 


- Kc 


(13) 
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-'c(  it  - J:(  

( 

The  coefficients  I and  F cf  ecuation  (12)  are  identical  to  £ and  1, 
respectively , except  v is  replaced  hy  v-j  in  the  arrunents  cf  the 
tessel  fnncticns. 


e Ticss  rate  or  ri.ow  at 


Series  Ir.pedance  cf  the  Annular  line.  7h 
the  reference  pr-essure  ?*  is  given  by 


= I (cu)  2rrdr 


ar.d  by  analogy;  tc  electrical  transTT.ission  line  theory 


(le) 

where  Z is  the  frequency  dependent  series  ir.pedance  per  unit  length , 

IF  is  analogous  tc  voltage  ar.d  C-  is  analogcus  to  c’urrent . 

Solving  for  (tu)  froc.  equation  (11),  perfom.ing  the  integration  cf 
equation  (1£),  and  using  equation  (16),  an  extression  for  Z is  obtained: 


„ iuc"  , 2ry. 


fF 


vr.ere 


where  Z is  shown  to  be  frequency  dependent. 


Shur.t  Adir.ittance  of  the  Annular  Line.  Fron.  electrical  trar.sTTiissior. 
line  theory. 


= YIF 


(20) 


where  j is  the  shur.t  aat.ittance  per  unit  length.  Upon  differentiating 
equation  (15)  wit;,  respect  tc  z,  using  equaticn  (1),  and  substituting 
into  equation  (20),  Y becoDes 


Fo 


Y = 


' r; 


00 

it 


r or 


(21) 
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t is  sr.cvr.  ir.  An  er-cix  i tr.a: 


££  U 

Y 


+ Ji 


5P 


/ C_  0 " 0 ■ 


(22) 


introducing  the  complex  time  notation  fcr  Sr/ot  and  ?T/;t  and  using  the 


relation 


cn  - 


ciT 


c e 


(23) 


equation  (22)  becomes 


^ jclF 

ot  V 


(•y-1)  e c-  AT 
1 + (y-1) *■ 


(24) 


where  £_  = 1/7  fcr  a perfect  gas.  Substituting  for  p Cp  AT/AP  from 
equaticn  (12),  equation  (24)  becomes 


?t  Y 


1 

/ 

1 - (y-1) 

EJo(  ^ 

u:  .^42 


FKo  ^ 


(25) 


Substituting  equation  (25)  into  equation  (21)  and  pei forming  the 
indicated  integration  yields 


y = 


1 + 


2it(1-y)  N 


T -* 


(2e) 
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ideal  gas,  = .-.-yv?*.  Tabulated  Bessel  functior.s  cf  interest  are 
giver,  in  [5,6]. 

Zcni valent  Zlectrical  Circuit . The  series  impedance  and  shunt 
admittance  may  be  expressed  as 

Z = R + jud  (26) 

Y = G - ja>C  (25) 

vhere  R,L,G,  and  C are  the  circuit  parar.eters  cf  resistance,  inductance, 
conductance,  and  capacitance,  per  unit  length,  respectively.  The 
equivalent  circuit  for  a transmission  line  of  dz  length  is  shown  in 
Figure  2.  The  circuit  parameters  are  all  real  numbers  and  are  dependent 
on  the  signal  frequencp-,  the  cross-sectional  area,  and  the  radii  of  the 
a.nnular  line  through  the  expressions  for  Z and  Y,  equations  (15)  and  (26). 

Characteristic  Imnedance  and  Propagation  Operator.  Zquations  (16) 
and  (20)  are  used  to  obtain  the  governing  differential  equations  for 
pressure  and  flow,  differentiating  equation  (16)  with  respect  to  z and 
using  equation  (2C)  gives 

ZY  (IP)  (30) 

82^ 
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he  renerel  soluticr.  fcr  the  t-ress’^e  is 


L?  - ^Z^  e ‘ “ -r  Z'-.e 


wr.ere  * , the  trcpagatior.  operator,  is  terineo  2S 


. =»'„;  = a - 


where  o is  the  attenuaticr.  per  ur.it  length  and  6 is  the  phase  angle 
i.'.crease  per  unit  length.  The  solution  for  the  vcluir.etric  flow  is 


- ■ ir 


wnere  Z;,  tne  oharaoteristic  ir.pedance,  is  defined  as 


Zo  = = 1/Yo 


and  Y;  is  the  characteristic  admittance. 

The  characteristic  impedance  and  propagation  operator  are 
generally  complex  quantities  and  are  convenient  for  describing  the 
pressure  and  flow  in  the  transmission  line.  An  excellent  review  of 
the  significance  and  use  of  Zq  and  T in  various  representations  of 


f ■ uid  lines  for  system  analysis  is  given  by  Goodson  and  Leonard  [7]. 


The  phase  velocity  is 


c = At  = 


2r.  6 


where  > is  the  wavelength  and  £ increases  2-n  radians  per  wavelength. 


The  phase  velocity  c is  in  general  not  equal  to  c^,  the  adiabatic 
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The  prcpagetior:  velocity  ratio  cr  waveler-ttr.  ratio  is 


c X 


The  oeper  attenuaticr.  in  distance  is 


a c=/u)  - a 


Tc  convert  the  attenuation  fror  nepers  to  the  ir-cre  corinor.  decibel, 
the  fcllovrint  ecuations  are  used: 


(aXj.)  db  = (logiQe)  a c^/o 


= 54.575  a Cg/ic  db  per 


(aX)  db  = 54.575  a/£  db  oer  X 


Characteristic  Frequencies.  Following  Nichols  [3]  it  is  convenient 
tc  define  two  characteristic  frequencies  which  will  aid  in  the 
presentation  of  the  results  and  facilitate  compariscn  of  the  annular 
line  results  with  those  cbtain=>G  for  circular  and  rectangular  lines: 


ro'd-r----^) 


5 r 


6rv>T  uj,, 

_ A _ v 

“T  = — %2 


where  r»'  = r^/rj,  is  defined  as  the  annular  radius  ratio.  The  viscous 
characteristic  frequency  is  used  to  fonr.  the  dimensicnless  frequency 
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The  ir.fluer*ce  cf  t:.e  fre.uency  ar.d  radius  ratic  cr.  ar.r.uiar  line 
perfornar.se  is  easily  seer,  in  Figures  3 thrcugh  6.  The  r^sculus  cf  the 
adr.ittarice  ratio  given  by  ecyuaticn  f35)  is  plotted  in  Figure  3 and  the 
phase  ar.gle  of  the  characteristic  admittance  is  plotted  in  Fig’ure 
It  car.  be  seen  that  at  high  frecyuer.cies  the  r.agnituce  cf  the  or.aracter- 
istic  ad.T.ittar.ce  approaches  that  cf  the  lossless  line,  vhile  tne  phase 
angle  approaches  cere.  At  low  frequencies  the  phase  angle  approaches 
degrees.  The  propagation  velocity  ratio  c/cg,  cr  wavelength  ratic, 
X/).g,  equation  (^2),  is  plotted  in  Fig-ure  3.  Fcr  a given  flow  area  arid 
frequency,  tne  propagaticn  velocity  decreases  with  increasing  radius 
ratio.  The  db  attenuation  per  line  wavelength,  equation  (‘•5),  is 
plotted  in  Figure  6.  For  given  fluid  properties,  flow  area,  and 
frequency,  the  attenuation  per  line  wavelength  increases  with  increasing 

rdGlUS  rdTlO* 

The  circuit  parameters  R,G,L,  and  C are  norm.alized  with  respect  to 
suitable  reference  values  and  the  results  plotted  in  Figures  7 through 


The  ratic  .R/Kyr.  is  plotted  in  Figure  7,  where  F-vr.  is  the  viscous 
laminar  steady  flow  (DC)  resistance  cf  the  a.nnular  line,  and  is 
defined  in  terms  cf  r*  [6]. 


8ru  l-t-r'-''^  1 


Equations  for  deterrr.ining  the  IC  properties  cf  circular,  annular,  and 
rectangular  lines  are  derived  in  Appendix  A.  As  the  frequency  u-  approaches 
zerc,  R approaches  Rv^i  however,  Ryr,  ^s  a good  approximation  for  R when 
ui/iuy  < 1 fcr  all  r”.  For  the  larger  values  of  r”  shown  in  Figure  7, 
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tr.e  l ir.ge  cf  ipilica.  ility  of  ccr.itar.t  r'-: Isoinct  r.:fel  iz  txtencef 
tc  higf.sr  values  of  u.  '^v 

riEure  5 shows  the  frecuer.cv  derer.dence  of  L-'l-  where  li  is  the 


adiaratic  ir.ductance.  The  tendency  cf  the  frecyaency  dependent  induct- 


ance to  artroach  its  acidbatic  value  at  high  frec.uenciec 
For  low  frequency  ratios,  say  u/u--.^  < 1,  the  inductance  is 
dependent  on  r*  and  essentially  independent  cf  ui/u..^.  The 
at  low  frequency  ratios  varies  between  apprcxir.ately  1.2 
depending  on  r*. 


is  clear 
prit.ari 
value  0 
and  1.33 


ly 

£ 


L/ 
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The  ratios  G/u,.-  Ca  arid  C/Ca  are  plotted  in  Figures  9 and  10.  In 
Fig'ure  9 the  conductance  is  shown  to  be  essentially  zero  when  < 1 

for  dll  r*.  For  values  cf  > 10  the  conductance  increases 

sigr.ificantly . This  increase  is  Ttore  pronounced  at  large  values  cf  r*. 
Figure  10  shews  the  general  thenriodynat.ic  behavior  cf  the  tiodel: 


C _ _ 1 (49) 

yF 

where  y is  the  adiabatic  exponent  (approximately  1.4  for  air  at  standard 
conditions)  and  n is  the  polytropic  exponent.  At  low  frequencies 
v/n  = 1.4;  thus,  n * 1,  which  is  the  isotheraal  exponent.  At  high 
frequencies  y/n  = 1.0,  thus,  n “ y»  the  adiabatic  exponent.  The  process 
is  therefore  isotherraal  at  low  frequencies  and  adiabatic  at  high 
frequencies.  As  r*  increases,  the  transition  from  isothermal  to 
adiabatic  behavior  occurs  at  higher  values  of 
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rif’,ur’<''  10.  Rntio  <A  F'rf'i^uency  Oepoinlf’iit.  to  A'l  i.nl'.it  i r C.ipnc  i t anco 


*r.  s'jr.r.ary , at  _cw  rretuenctea  tr.e  et/jiva^s 


circuit,  see 


Figure  2,  ccnsists  essentially  of  constant  value  elements  of  resistance, 
inductance,  and  capacitance  and  no  conductance.  At  high  frecuencies , 
the  inductance  and  capacitance  apprcach  their  adiabatic  values;  hcv;ever . 


cn  trecuencv  and  radius  ratio. 


t^rr.tariscn  with  Circular  and  Fectancular  lines 


The  results  obtained  fcr  the  annular  line  are  now  ccr.pared  with 
these  obtained  fcr  circular  and  rectancular  lines.  The  cemoarison  is 
r.ade  cn  the  basis  that  the  respective  crcss-secticnal  areas  are  equal. 

Generally,  the  annular  results  approach  those  fcr  circular  lines 
as  the  inner  radius  cf  the  annulus  approaches  zero  and  the  losses 
attributable  to  the  inner  wall  becoTne  small.  It  can  be  shown  that  the 
circular  line  velocity'  and  temperature  profiles  of  .Nichols  [3]  may  be 
obtained  directly  from  equations  (11)  and  (12)  for  the  special  case  cf 
the  annular  line  with  r^  = 0.  The  Bessel  function 


-IF..-) 


is  infinite  at  the  origin  and  ceases  to  be  a solution  in  equations  (11) 
and  (12).  Thus,  when  r^^  = 0,  the  constants  D and  ? are  identically 
zero.  Applying  the  bcundary  conditions  of  equation  (1C)  at  the  outer 
wall  r = Tq  and  obtaining  new  relationships  for  B and  E leads  to 
Kichcls ' expressions  for  the  velocity  and  temperature  profiles.  In 
all  of  the  figures  the  curves  for  r'“’  = 0 are  identical  with  those  of 
Kichols  for  circular  lines. 
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The  anr.'jlar  resu, 
lir.es  as  r*  1.  The 
ratio  for  the  ar.nular 
d65ir.€d  hsr's  3S  AR>*  - 


.ts  converge  to  those  of  Schaedel  [*•]  for  rectangular 
comparison  is  made  by  defining  an  equivalent  aspect 
line.  The  aspect  ratio  of  the  rectangular  line  is 
b/h  > 1.  For  equal  rectangular  and  annular  flow 


areas , 


bh  = r(ro^-ri^)  = -(r^+r^ ) (rj.-r4  ) 


(50) 


where  choosing  b = rCr^-^rq)  and  h = r^-rj^,  the  annular  aspect  ratio  is 
defined  as 


. rCltr^O 
(rl-ri)  (l-r--) 


(51) 


Since  the  annular  line  side  walls  cf  dimension  h dc  not  physically 
exist  and  are  net  considered  in  the  theory,  greater  losses  are  expected 
in  rectangular  lines  when  compared  to  annular  lines  of  the  same  area 
arid  aspect  ratio.  However,  as  aspect  ratio  becomes  large,  the  losses 
attributable  to  the  short  wells  of  the  rectangular  line  become  small 
compared  to  the  losses  on  the  long  walls.  Consequently,  the  results  for 
annular  and  rectangular'  lines  converge  wioh  increasing  aspect  ratio. 
Figure  11  is  a coriparisen  of  the  attenuatior.  per  unit  wavelength  (aX) 
cf  annular  and  rectangular  lines  where  here  is  defined  to  correspond 
tc  tnat  used  by  Schaedel  [»•];  i.e.,  Wy  = ^v/A.  Cf  interest  is  the 
rather  close  agreement  at  moderately  low  aspect  ratios.  Agreement 
within  10%  is  achieved  for  all  frequency  dependent  properties  when 
A?.  > 7 (r*  > 0.38).  Thus,  either  model  will  give  approximate  predic  ions 
of  the  perfcrm.ar.ee  of  the  other  model  when  AR  > 7. 
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The  differences  in  the  circular,  rectangular,  and  atinular  line 
results  fcr  equal  areas  are  due  ir;ainly  to  differences  in  the  steady- 
flow  (11)  resistances.  The  differences  in  the  TC  resistances  are 
illustrated  by  the  DC  resistance  ratios  plotted  in  Figure  12.  The 
viscous  resistance  per  unit  Tength  of  circular  lines  is  the  classical 
.Hagen- Peis euille  flow  resistance 


where  Rvj-.  is  given  by  equation  (46).  The  ratios  Rvr'’^'Vc 
are  alsc  given,  where  is  the  viscous  resistance  cf  rectangular  lines 
given  by  Schaedel  [4].  The  ratios  Rvn/^vc  Rvr’^^-vc  nearly 
linear  functions  of  AP,  fcr  AR  The  convergence  of  to  unity 

for  large  AR  is  easily  seen. 

In  s'jtjriary,  the  sTcail  sigTicl  perfomance  of  annular  lines  has  been 
presented  arjd  includes  the  perfort.ance  of  circular'  lines  as  a special 
case  (r*  =0).  The  small  signal  perfcm.ance  of  annular  arid  rectangular 
lines  converge  with  increasing  aspect  or  radius  ratio.  Useful  approxi- 
mations of  rectangular  line  performance  ii.ay  be  made  using  the  annular 
line  results  at  higher  aspect  ratios,  or  vice  versa. 

It  is  important  to  note  that  whereas  the  namerical  results  are 
presented  for  air  at  80^F,  any  other  fluid  could  have  been  used.  If 
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t.'e  fluid  is  considered  incor.pressibie , tde  results  are  easily  adapted, 
and  considerably  simplified  by  assuming  v = 1,  and  by  taV.ing  the  speed 
cf  scund  as  c = where  Yi  is  the  bulh  ncdulus  of  the  fluid.  It  is 

noted  that  the  s.fur.t  admittance,  equation  (2i),  becomes  a p-ure 
capacitance  of  constant  value  when  y = 1.  Tne  limitaticns  given  by 
Iberall  [1]  apply  tc  the  annular  line  also; 


rr-r.-  > — (=  1.6*1C“^  inches,  for  air) 

Ca 


u,  < (=  7.7x10-  rad/sec,  for  air) 

V 


The  annular  line  model,  as  well  as  the  circular  and  rectangular  models 
of  [3]  and  [u],  respectively,  are  restricted  to  lam.inar  flow. 
Discrepancies  are  to  be  expected  for  turbulent  mean  flow  conditions. 
The  case  of  turbulent  miean  flow  is  considered  in  Section  V. 


High  and  Low  Frecuenoy  Approxim.ations 

Useful  high  and  low  frequency  approxim.ations  cf  the  system,  variables 
are  most  readily  obtained  from,  arialysis  of  the  equivalent  circuit 
components  in  conjunction  with  the  applicable  Bessel  function 
apprcxim.ations . 

For  the  low  frequency  approxim:ations  <<  1),  the  equivalent 

circuit  consists  of  the  constant  value  elements  Rv^,  and  , with 

= C.  Here  Ci  is  the  isotherm, al  capacitance,  = YCg.  The  low 
frequency  series  im.pedance  and  shunt  admittance  are  respectively. 


+ jwLvjj 


'^-n  = + juCi  » jwYCa 
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using  the  inequality  u.Lv^/?-Vp_ 

Ic/W  fretuencv  tharacteristic  adr. 


1 and  equations  ' I"  ) arjd  (5-',  the 


and 


cor.' 


YCau  ^ 
irvp 


(Itj) 


(55) 


r.  z 


-Ivy'll 


a + ji 


•''VpCgu 


(1+:) 


(Si) 


It  is  seen  that  the  phase  lag  S is  approxirately  nurerically  equal  to 
the  attenuation  a,  and  that  the  phase  angle  of  the  adraittance  is 
approxirately  r/4  radians.  For  < C.l,  equations  (55)  and  (55) 

give  " and  Yj  within  one  percent  of  the  exact  value,  for  all  r*'. 
Expressions  for  Fvn  =“d  Ly-  are  derived  in  Appendix  A. 

The  high  frequency  approximations  are  derived  in  Appendix  E using 
the  standard  asymptotic  expansions  of  the  Bessel  functions  [5,5], 
along  with  other  approximations.  For  w/uy  >>  1,  the  inequalities, 

R/iuL  <<  1 and  3/u;C  <<  1,  apply.  The  high  frequency  series  impedance 
and  shunt  adm.ittance  become 


Z; 


3ui-ar. 


! 


1 + 


rrc 

^n 


*1  '[  (r-*-l)  + :’(r'=+l)] 


± 


(57) 


•Hn  = 3u>ar. 


[(r>'!-i)  + 


j(r»''+l)] 


(56) 


Tl-,e  high  fr-.equency  propagation  operatcr  and  characteristic  adr.ittance 
become 
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'h  = - :£h  = ^ ^ 


l+r*' 

l-rs'! 


^ ww . + 

,J 


^55) 


1 -’u’ 

^ n 

• a 

1 * -iL \ 

“ \ (i.  L u C / 

' 1 

(£C) 


It  is  cf  interest  that  the  extressicn  fcr  cii.;  in  equation  (55)  for 
annular  lines  is  identical  to  that  citained  ty  Karan  and  Fra-nke  [5]  for 
circular  lines  for  the  snecial  case  r*  = 0.  The  extressicn  for  £v  is 
identical  to  that  given  in  [5].  It  should  therefore  be  oossible  to 
obtain  a solution  in  the  time  domain  for  annular  lines  using  Karam's 
approach  [15, 2C].  The  frequency  to  time  domain  transformations  fcr 
circular,  rectangular  and  annular  lines  are  covered  in  Section  III. 


Ixperimental  Confirmation  of  Annular  Theory 

The  experimental  verification  of  the  annular  line  theory  was  the 
primary  purpose  of  a K.S.  thesis  completed  by  Captain  G.  R.  Farney 
at  the  Air  Force  Institute  of  Technology  [10].  Farney  measured  the 
small  signal  amplitude  and  phase  shift  frequency  response  of  a blocked 
annular  pneumatic  transm.issior.  line  at  five  different  radius  ratio.. 

The  outer  line  was  approximately  three  feet  long  and  had  an  internal 
dic.meter  of  0.167  inches.  The  five  radius  ratios  were  achieved  by 
inserting  and  centering  five  different  sized  rods  within  the  basic 
line.  Farney  observed  that  the  experimental  gains  were  within  ± 1 db 
and  the  phase  shift  angles  within  il5  degrees  of  those  predicted  by 
the  annular  line  theory.  Further,  the  predicted  analytical  results  fell 
within  the  measurement  capability  of  the  experimental  instrumentation  used . 
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A ^6£cnp“iwr.  C-  *r.^  sxp6rir.tr.Pal  apparsi'is  arc  tssP  rrX'Csd'irss  ussc 


is  giver,  ir.  [iC]. 


ec’-atior.s  are  required  tc  acccur.t  fcr  t'r.e  Eer.dirg  ar.d 


receiving  cor.citicns . These  are  stand; 


ecuations  or  eiectrica- 


tra.nsr.ission  air.e  theorp',  availalle  ir.  r/a.ty  references  such  as  [11] 
ar.d  are  repeated  herein  for  the  convenience  of  the  reader. 

For  a giver,  line,  the  transfer  pressure  gain  is  [12]: 


g = 'r^l 


a; — cash  r i 

■^S 


^r  . . 

TTo-  Sinn  i Jc 

4-0 


where  r and  s denote  sending  and  receiving  locations  and  1 is  t.ne 
distance  between  these  locations.  The  sendir.g  ir.pedar.ie  is 


Zj-.  cosh  n + Z;  sinh  Ti 
Zq  cosh  n + Z.-.  sinh  Ti 


For  a blocked  line,  the  receiving  impedance  Zj-  is  infinite,  thus 


The  phase  lag  angle  between  sender  ar.d  receiver  is 


fcl  = arctan 


U) 

Real  (g) 


The  pressure  gain,  equation  (£1)  converted  to  decibels,  and  the  phase 
lag  ar.gl.e,  equation  (6-+)  in  degrees,  are  plotted  as  functions  of 
excitation  frequency  in  Figures  13,  1*^,  and  15,  fcr  annular  radius 
ratios  cf  0.C7,  3.23,  ar.d  0.6*’,  respectively.  The  experimental  data 
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rigure  1 


Correlation  of  Experimental  Fesults  with 
Annular  Line  Theorv 
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Figure  15, 


Correlation  of  Experimental  Pesults  with 
Annular  Line  Theory 
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Steady 

flow,  it 

has  been  particular!; 

ices  cf 

lines  of 

ncr.circular  cress  se 

useful  to  correlate  tr.e 
tiers  with  the  resistarce 


of  circular'  lines  through  the  classical  hydraulic  diareter.  Sir.ilarly, 
in  oscillatory  flow,  it  would  appear  useful  tc  he  alle  to  correlate  the 
frequercy  response  of  ror.circular  lines  to  the  response  of  circular 
lines.  Besides  the  corr.pariscns  r.ade  herween  rectang'ular  and  circular 
lines  [4]  and  between  rectangular,  annular  and  circular  lines  herein, 
ether  attempts  have  been  r.ade  to  predict  the  response  of  rectangular 
lines  using  the  concept  of  the  hydraulic  diameter  [13,14], 

In  this  section,  a dimensionless  frequency  is  defined  which  leads 
tc  a generalized  presentation  of  the  frequency  dependent  response  of 
lines  cf  circular  and  noncircular-  cross  section.  The  dimensionless 


frequency  in  turn  leads  tc  a characteristic  radius  for  the  noncircular 
line  which  is  rcughly  analogous  to  the  classical  hydraulic  radius  used 
in  steady  developed  flow  for  noncircular  lines.  It  is  shown  that  the 
fr-equency  dependence  cf  noncircular  lines  may  be  very  closely  approxi- 
mated when  the  characteristic  radius  is  used  together  with  the  circular 
line  theory  cf  i’ichcls  [3],  The  results  are  compared  with  those 
obtained  by  .Healey  arid  Carlson  [13]  for  rectangular  lines.  In  [13],  the 
classical  hydraulic  radius  together  with  Nichols'  theory'  was  used  to 
predict  the  response  of  rectangular  lines.  The  response  predicted 
using  the  characteristic  radius  defined  herein  is  shown  to  be  clearly 
superior  to  the  response  predicted  using  the  hydraulic  radius.  Results 
are  presented  for  both  rectangular  and  annular  lines  over  a wide  range 


cf  dipect  and  radius  ratio,  respectively.  In  addition,  ar.  approach  tc 
iTTitedance  r.at chirr  based  or.  ecuivalent  circular  lines  is  t resented. 


Frequency  .Response  Solutions 

Several  dimensionless  frequencies  (cr  frequency  ratios)  have  beer, 
defined  in  order  tc  present  generalized  frequency  response  solutions 
for  circular  lines,  l.'ichols  [3]  showed  that  the  frequency  response  cf 
a circular  line  is  a functicr.  of  the  dimensionless  frequency  ratio 

where  u is  the  angular  excitation  frequency  and  is  a viscous 


cnaractercs' 


cequencj'  aefinec  as 


‘“V  = i: 


.Ry  is  the  laminar  DC  viscous  resistance  per  unit  line  length  and  is 
the  adiabatic  inertance  per  unit  line  length.  For  any  line,  and 
are  defined  as 


^ ft] 


For  a circular  line,  Ry  and  Wy  become 


'c  , 2 


‘"^'c  = T 
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Other  forr.ulatioris  cf  the  ciKensior.less  frequency  have  also  heer, 
defir.ed  for  the  circular  line;  for  example , BroKn  et  al.  [15]  used 
r.  = L.r^/v  = 8u/a.'vj.,  while  Golds chir.iec  [15]  used  Sh'  = dur^/v  = 32 
In  ail  of  these  cases,  however,  the  dimensionless  frequencies  are 
proportional  cc  cA/v. 

Schaedel  showed  that  the  response  of  recr angular  lines  of  the 
same  aspect  ratio  is  a function  of  aiA/4v,  and  it  has  teen  shown  herein 
that  the  response  cf  concentric  annular  lines  cf  the  same  radius  ratio 
is  a function  cf  loA/Btv.  Tp'pical  results  for  the  attenuation  per  line 
wavelength  of  rectangular  arjd  annular  lines  are  plotted  in  Figures  16 
and  17,  respectively.  Circular  line  results  are  included  in  Figure  17 
as  a special  case  of  the  annular  line  with  the  inside  radius  r^  = 0. 

The  dependence  on  aspect  ratic  and  radius  ratio  as  well  as  signal 
frequency  is  clearly  shown. 

Generalized  Response  of  Koncircular  Lines 

The  aX  curves  shown  in  Figures  16  and  17  for  different  aspect 
ratios  and  radius  ratios  all  have  the  sa.me  general  shape.  Thus,  it 
would  seem  possihle  that  the  results  for  different  aspect  or  radius 
ratios  could  he  generalized  and  related  tc  a single  curve.  This  can  be 
accomplished  approximately  by  using  the  dimensionless  frequency  defined 
by  equation  (65)  with  Ry  taken  as  the  resistance  cf  the  particular 
cross  section  being  considered. 

A better  approach  is  to  include  the  inertance  effects  in  the 
dimensionless  frequency  as  well  as  the  resistance  effects  and  to  present 
the  results  as  a function  of  the  dimensionless  frequency  where  01^ 

is  the  radiaij  frequency  defined  by  Kichols  [3]  as 
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rigurp  J.7.  Alleiiuatitm  Ppr  l.irip  Wavelptigtli  for  AniiuJrH-  l.iiirr 


’^c  * 

-V 
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Here  Hv  is  defined  by  ecyjation  (£?\  and  !•/,  tr.e  lcv;-frev-^" "v  '^- 


:us  ir.ertance,  is  defined  by 

-■  = [j ! (t/  ^ ^ 


where  depends  on  the  ianinar  fully- developed  velocity  distritution. 

The  ratio  l./u)c  includes  both  viscous  and  inertial  effects  at  low 
frequency.  Consequently,  this  ratio  ir.cre  accurately  represents  the 
fluid  line  at  lew  frequency  ratios,  and  gives  analytical  results  which 
are  in  closer  agreerr.ent  with  both  the  theoretical  and  excerinental 


.ar  and 

noncircular 

lines . 

n Figur 

'€S  15  and  17 

are  retlcttec 

Ir. 

The  families 

of  curves  for 

di: 

aspect  and  radius  ratios  are  now  reduced  approxir.ateiy  to  a single 
curve,  ry  incorporating  the  geometric  dependence  cf  aspect  and  radius 
ratio  in  lines  cf  different  sizes  and  cross  sections  are  shown  tc 
have  approximately  the  same  generalized  frequency  tehavior.  The 
circular  line  attenuation  is  also  included  in  Figure  15,  where  for 
the  circular  line  is  the  same  as  that  given  by  Nichols  [3];  i.e.. 


C '"•c‘  61IV  3 


- -ac 


= 4 “Vc 


Figure  15  shows  aX,  the  attenuation  per  line  wavelength.  In  order  to 
obtain  a,  the  attenuation  per  unit  length,  it  is  necessar\’  to  know  X 
the  wavelength  in  the  line  for  each  ratio  iii/w(|,.  Curves  of  X/Xg  vs 
u)A/4v  for  rectangular  lines  are  given  in  [4],  and  are  not  reproduced 
herein.  However,  for  annular  lines  curves  of  X/X.  vs  uiA/Bttv  are  given 
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Nur.er  ical  Results  and  lis  cuss  ion 

The  results  in  Fig'ures  1£  through  1?  were  ccr.puted  assuming 
standard  properties  of  air  at  SC  F (Table  I).  The  numerical  differences 
in  in  Figure  IS  for  the  range  of  aspect  or  radius  ratios  shovrn  in 
Figures  15  and  17  were  small  and  are  not  plotted  in  Figure  IS.  The 
maximum  difference  in  a>.  between  circular  and  rectangular  lines  was 
1.05  db  at  u/u:^  = 6.C  for  a = 10.  The  miaximum  difference  between 
annular  and  circular  lines  was  i.4C  db  at  = 5.0  for  r"  = O.B.  For 

u 7he  difference  was  -ess  than  0.0?  db,  and  for  w/wc  > 10  the 
difference  was  less  than  1.25  db. 

The  maximu.m  difference  in  between  circular  and  rectangular 

lines  was  C.03E  at  u'/uc  = 2.0  and  a = 10.  The  maximum,  difference 
between  circular  and  annular  lines  was  0.046  at  w/wc  = 3.0  and  r"  = 0.6. 
For  < 0.1  the  miaximum.  difference  for  bcth  rectangular  and  annular 

lines  was  less  than  0.01"^  and  for  u/ui^  > 10  the  maximum,  difference  was 
less  than  0.025.  The  m.aximum.  departure  cf  fromi  that  predicted  for 

t.ne  circular  line  for  the  range  of  aspect  and  radius  ratios  considered 
occurred  for  the  annular  line  of  radius  ratio  r*  = 0.8,  as  shown  in 
Figure  19.  All  other  configurations  are  nearer  tc  the  circular  line 
results . 
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It  is  noted  that  whereas  a)  vs  w/u,  gives  the  tightest  rroupir.z  of 
the  aV  curves,  \/>.^  vs  gives  a sligr.tly  better  grouting  cf  the 

X/^3  curves,  where  is  defined  by  eoyuatior.  (65). 


Setero-ination  of  the  Characteristic  Frequency 

The  deterr.inaticr.  of  values  of  for  ncncircular  lines  reduces 
essentially  to  obtaining  values  of  Ry  and  by.  Analytical  expressions 
cf  r-v  and  Ly  -or  circular,  rectangular'  and  annular  lines  are  developed 
arid  discussed  in  Appendix  A.  A convenient  rr.ethod  for  obtaining  F.y  is 
tc  relate  the  of  the  line  with  the  ncncircular'  cross  section  to  that 


^^V  _ y 

•^Vc  A^ 


where  Kj;_  for  rectangular  lines  and  annular  lines  depends  on  aspect 
ratio  and  radius  ratic,  respectively.  Values  of  .Kp_  have  been  calculated 
[**,17]  and  a few  are  giver,  in  Table  III.  The  viscous  inertance  by  is 
obtained  fron.  values  of  Kl,  equation  (71),  for  each  geometric 
configuration.  Values  of  have  also  been  calculated  [‘*,17,15]  and  a 
few  are  given  in  Table  III. 

Starting  with  equation  (75),  a convenient  equation  for  computing 
is  easily  derived  as 


^ _8itv 

Kl  a 


For  circular  lines:  = 1 and  KLj,  = 4/3;  thus  = Girv/A,;.  Values 

of  for  noncircular  lines  can  easily  be  obtained  using  values  from 
Table  III  or  using  the  equations  or  data  from  the  references  [4,17,18]. 


Table 


The  steady  flow  iat.ir.ar  resistance  Rv  does  not  depend  on  frequency; 
consequently,  Ry  ir.ay  be  expressed  in  terir.s  of  the  classical  steady- 
flow  parameters: 


Rv  = fRe 


IL-h^A 


(75) 


where  Iv^  is  the  conventional  hydraulic  diairjeter,  f is  the  friction 
factor  and  Re  is  the  Reynolds  number.  For  a circular  line,  Dy,  = Ec 
and  fRe  = 6*^;  therefore, 


■'Vr 


= e*4 


(76) 
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lividirig  equation  (75)  ty  equation  (7£) 


(77) 


where  fne,  iv,  and  A .are  determined  for  the  noncircular  cross  section 
being  considered.  Selected  values  of  f.Re  are  given  in  Tao.le  III  for 
rectarigular  and  annular  cross  sections.  The  relaticnship  between  K.p 
and  f?.e  is  obtained  from  equations  (73)  and  (77) 


% 


fP.e 

6u 


A 


(76) 


where  Kq  is  a geometric  shape  factor  defined  as 


K 


T-  ^ ^ 

-'c  A _ 4A  _ ^ w 


(75) 


Kq  is  a constant  for  each  geometric  cross  section.  It  is  important  to 
note  that  when  Dj-  = the  areas  of  the  noncircular  cross  section  A and 
equivalent  dia-meter  circle  Ac  ate  not  equal  and  rq:  = k/f-.Q.  Likewise, 
when  A = the  hydraulic  diamieters  are  not  equal  and  Kq  = Bc^/Dv^. 
Expressions  for  Kq  are  given  in  Table  IV  for  circular,  rectangular  and 
annular  cross  sections.  Selected  values  of  are  included  in  Table  III 
fcr  rectangular  and  annular  cross  sections. 
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TA2LE  IV 

leoinetric  Cross  SecticTi  Factor 


Use  cf  the  Hvcraulic  Diar.eter 

Results  hasec  cr,  hydraulic  diat.eter  car.  he  used  to  tredict  line 
perfcrTiance  for  STTiall  signal  oscillator^’  flo’A-.  At  a giver,  signal 
frequency,  Figure  1£  shows  that  for  the  sane  attenuation  on  lines  of 
different  cross  section,  the  itust  he  the  sanie.  The  use  cf  equivalent 
hydraulic  dia-neter  for  a noncircular  line,  however,  leads  to  an  Uj-  that 
is,  in  general,  not  equal  to  This  will  now  be  illustrated  fcr  a 

rectar.gular-  line. 

Fror.  equation  (74),  the  ratio  u)cv./wcc  is 


i^Cr  4 A(- 

“C(-  3 Ay. 


(80) 


By  substituting  for  froir.  equation  (78),  the  ratio  becoraes 


Wf'  ^ Xi c,  - n^ 

“Cc  3 64  K-  r).2 

‘*r 

For  = I!^,  equation  (81)  reduces  to 


(81) 


I 
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3K-,. 


eu 


(S2) 


Typical  values  of  and  as  a function  of  aspect 

ratio  are  shown  in  Tahle  V using  the  values  of  fr^^r  £nd  from 
Talle  III. 


TAEhl  V 

Frecyaency  Ratio  Rarat.eters  for  Rectangular 
Lines  fcr  Dn^  = 


b/h 

uicr/wcc 

v'u-cp/wcc 

0. 86 

0.93 

2 

0.S6 

0.96 

3 

1.09 

1.  04 

4 

1.18 

1.09 

5 

1.25 

1.12 

10 

1.4c 

1.19 

The  ratio  car.  be  interpreted  approximately  as  the  ratio  of 

the  viscous  attenuation  cf  rectar.gular  to  circular  lines  of  the  sane 
hydraulic  dianeter  at  a given  frequency.  The  values  of  given 

in  Table  V confirm,  the  attenuation  results  of  Healey  ar^d  Carlson  [13] 
arid  explain  the  differences  they  found  between  the  rectangular  line 
solutions  and  circular  line  solutions  cf  the  same  hydraulic  diairieter. 
VTnen  using  circular  line  equations  to  predict  rectangular  line  results, 
the  differences  cazi  be  avoided  if  the  rectangular  line  friction  factor 
(fj.rex-)  and  inertance  Lv^,  are  included  in  the  calculations.  This  pcint 
will  be  made  clearer  in  the  next  few  paragraphs. 
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I/cter~iri£~icr.  cf  the  Characteristic  Radius 

It  is  irr.pcrtar.t  tc  note  that  characterizes  the  frequency 
dependent  line  perfcrr.ance , and  that  tc  coir.pute  u;,  only  the  hinet.atic 
viscositt',  line  geor.etry,  and  steady  flcv;  lat.inar  velocity  profile  need 
oe  known.  rresur.ir.g  that  u-  = is  known  for  a line  cf  arritrar-y  but 
constant  axial  cross  section  an  equivalent  circular  line  containing  the 
sa.7:e  fluid  is  sought.  The  characteristic  frequency  of  the  circular 

• ^ rife 


£v 


Setting  iuqq  - the  radius  of  an  equivalent  circular  line,  called 

the  characteristic  radius,  is  easily  found  to  be 


Since  no  additional  restrictions  have  been  placed  on  the  arbitrary  line 
(i.e.,  neither  the  flow  areas  nor  hydraulic  diameters  are  assuiTjed  equal), 
equation  (84)  implies  that  for  each  arbitrary'  line  there  exists  a unique 
equivalent  circular'  line.  The  equivalent  circular  line  w’ill  have 
approximately  the  same  normalized  frequency  dependent  behavior  as  the 
arbitrary  line  since  the  characteristic  frequencies  are  equal.  The 
characteristic  radius  is  roughly  analogous  to  the  hydraulic  radius  used 
in  steady  flow  problems. 


«< 

1 

5 


I 


-Crr.i,  ar 


of  '^narccteristic  ar.c  nvcrau-ic  .^acol 


Heaiey  ar.d  Carisors  [13]  used  the  nycraulic  radius  (hasec  on  tne 
iu*ic  diar..eter)  for  rectangular  lines, 


r-^ 


(83) 


together  with  Kichcls'  equations  for  circular  lines  [3]  to  estiir-ate  the 
frequency  response  of  blocked  pneun.atio  lines  of  rectang-ular  cress 
section.  These  results  were  then  compared  with  the  results  predicted 
by  the  rectangular  line  theory'  of  Schaedel  [^^3  and  with  their  own 
experimental  data.  The  experirriental  data  and  predictions  using 
Schaedel 's  theory  were  in  excellent  agreement  for  all  aspect  ratios 
tested  (a=l,2,3,b,5).  Using  r>-,  and  Nichols'  equations  for  circular 
lines  yielded  good  agreement  with  the  data  and  with  Schaedel 's  theory 
for  aspect  ratios  of  2,  3,  and  i-;  however,  at  a=l  and  a=5  significant 
differences  were  observed  between  the  predicted  downstream,  amplitude 
response  using  rj,  and  the  mieasured  response. 

Using  the  characteristic  radius,  equation  (6-),  together  with 
Nichols'  theory,  the  ncrTialized  attenuation  aX  vs  u was  obtained  for 
aspect  ratios  a = 1,  3,  5,  and  10  and  the  results  plotted  in  Figure  20. 
For  comparison,  predictions  using  rrj  and  Nichols'  theory  as  in  [13], 
and  using  Schaedel 's  theory  are  also  plotted.  The  differences 
between  Schaedel 's  results  and  those  using  r^  are  so  small  that  they 
cannot  be  seen  on  Figure  20.  It  may  be  concluded  that  the  use  of  r,; 
gives  much  better  agreement  with  Schaedel' s theory  and  consequently 
with  the  experimental  data  of  [13],  than  does  the  use  of  rq,.  It  can 
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radius  predicts  toe  r.uch  atter.uatiori 


aisc  te  seen  teat  the  use  of  hydraulic 
per  line  wavelength  fer  a=l  and  net  enough  for  ai 
increases  with  aspect  ratio  aheve  a=3. 

A sir.ilar  comparison  is  presented  for  annular 
0.1,  C.£,  ana  C.S  an  Tagure  2±.  At  r"'=0 , r^^  = r^ 


The  difference 


fer  circular  lines  is  recovered.  For  r->0,  the  attenuation  predicted 
using  r-  and  Nichols'  thecri'  is  less  than  that  predicted  using  either 
r-  and  Nichols'  theor\’  cr  the  exact  annular  thecr'y  given  in  Section  II 
while  the  latter  two  approaches  are  in  excellent  agreenient.  The 
difference  in  predicted  perforniance  based  or.  r^  and  r;..  continue  to 
diverge  with  increasing  radius  ratio. 

It  is  of  interest  to  ccm.pare  the  hydraulic  radius  rr.  and  the 
characteristic  radius  Vq  of  a rectarjgular  line,  for  a practical  range 
of  aspect  ratio,  Figure  22.  The  rwo  approaches  yield  the  sar.e 
equivalent  circular'  line  at  a=2.3  (ri;=r-,,),  and  are  nearly  equal  ever 
the  rarjge  2 < a s 4.  However,  the  divergence  of  r^  and  rv,  away  froir. 
a-2.3  is  clearly  seen.  The  circular  radius  required  to  give  the  sane 
flow  area,  denoted  r^t  is  defined  as 


= /A/t 


(5e) 


and  is  used  tc  ncrrrialize  both  rv  and  r^.  Similar  results  are  plotted 
for  annular  lines  in  Figure  23.  It  can  be  seen  that  the  ratio  rc/ry, 
decreases  rapidly  from  unity,  where  r»''=0,  and  = rj._  = Tf.  Th<=«  area 
of  the  equivalent  circular  line  using  the  characteristic  radius  is 
apprcximately  60%  of  that  obtained  using  the  hydraulic  radius  for 
> 0.2. 
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In  principle,  if  the  frequency  dependence  of  lines  of  arbitrary 


axial  cross  section  may  be  correlated  as  above  using  then  an 
approximate  impedance  match  may  be  fo'jnd  when  it  is  necessary'  to  join 
lines  of  different  axial  cross  sections,  as  discussed  in  the  next  few 
paragraphs . 

Impedance  Matching 

It  is  of  practical  importance  to  terminate  a transmission  line  in 
its  characteristic  impedance  or  admittance  since  such  a termination 
implies  no  reflections.  Minimizing  reflected  signal  strength  is 
synonymous  with  increasing  the  efficiency  of  the  line  in  transm-itting 
the  desired  signal.  Thus,  it  is  appropriate  to  discuss  impedance 
matching  when  lines  of  different  axial  cross  section  are  to  be  joined. 
Such  circum.stances  occur  quite  often  in  fluidic  applications  where 
circiilar  lines  feed  fluidic  components  whose  basic  flow  regions  are 
rectangular  due  to  laminate  Techniques  in  component  fabrication.  The 
joining  of  flow  regions  of  different  axial  cross  section  may  also  be 
found  in  aircraft  hydraulic  systems. 

If,  for  the  present,  we  ignore  the  junction  itself  and  concentrate 
on  matching  the  characteristic  impedance  of  transmission  lines  of 
different  cross  section,  it  wo'uld  seem  appropriate  to  require  the 
characteristic  frequencies  to  be  matched. 

The  modulus  and  phase  angle  of  the  characteristic  admittance  ratio 
for  circular,  rectangular  and  annular  lines  versus  oj/ai^  are  plotted  in 
Figures  24  and  25,  respectively.  Again  it  is  seen  that  lines  with  the 
same  exhibit  approximately  the  same  frequency  dependence.  The 
admittance  ratio  of  rectangular  lines  is  within  4%  of  that  of  circular 


Circular,  Rectangular  and  Annular  Lines 


lase  Anglp  of  Characteristic  Admittance  for  Circular,  Rectangular  and  Annular  Lines 
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lines  at  aspect  ratios  a < .10,  and  the  adtittance  ratio  of  annular 
lines  is  within  6%  of  that  of  circular  lines  for  radius  ratios  r*  < O.S. 
At  lower  aspect  or  radius  ratios,  these  percentages  decrease.  The 
phase  angle  of  the  adnsittance  vs  w/loq  for  both  rectangular  and  annular 
lines  is  within  1.5  degrees  of  that  of  circular  lines  for  the  sane 
ranges  of  aspect  and  radius  ratios. 

It  is  important  to  note  that  the  modulus  of  the  characteristic 
adiiiittance  has  been  nonrialized  with  respect  to  that  of  a lossless  line 
o.f  the  same  cress  section.  In  general  terms 


lYo  1 

= f(u/(t'  ) 

I^SOl 

or 

lYol  = lY3ol-f(i./a)c) 

(87) 

and  equations  (87)  are  valid  for  all  lines.  It  follows  that  if  two 
lines  have  the  sane  characteristic  frequency,  then  for  a given  excita- 
tion frequency  u the  lines  will  have  approximately  the  samje  admittance 
ratic  but  not  necessarily  the  same  admittance.  From  the  definition  of 
^so  ’ 


^so  = ^ 


^ P*ca 


(36) 


and  YgQ  is  obviously  area  dependent.  Therefore,  setting  the  character- 
istic frequencies  equal  does  not  insure  that  the  characteristic 
admuttances  are  equal.  If  the  conditions  that  Yq  and  must  be  the 
same  are  imposed  on  two  lines  carrying  the  same  fluid,  then  the  flow 
areas  must  also  be  the  same.  These  conditions  are  in  general  incompat- 
ible except  for  the  special  case  when  the  two  lines  are  identical.  The 


66 


I 


conditions  ajno'jnt  to  requiring  that  = r^  when  lines  of  oircular  and 
noncirc'olar  cross  seotion  are  to  be  joined.  Fror.  Figure  23,  it  is  seen 
that  for  annular  lines  only  when  r*  = C;  for  rectangular  lines, 

Figure  22,  this  condition  is  never  achieved.  It  is  concluded  that  the 
characteristic  admittance  (or  impedance)  of  lines  of  different  cross 
section  carrying  the  same  fluid  cannot  be  matched  at  each  and  every 
physical  frequency  w.  The  closest  match  to  a circular  line  is  obtained 
when  r6=rc.  (r*‘  ->•  0 or  a 1,  for  annular  or  rectangular  lines.) 

In  practice,  the  most  common  constraint  encountered  is  the 
specification  of  a steady  mean  flow  through  which  the  signals  are  to 
be  transmitted.  This  constraint  is  equivalent  to  requiring  that  the 
flow  continuity  is  maintained.  It  follows  that  for  a given  fluid  two 
lines  of  the  same  flow  area  have  the  same  lossless  admittance.  At  a 
specified  frequency  ratio  w/wc,  tlie  two  lines  will  also  have  approxi- 
mately the  same  frequency  dependent  characteristic  admittance.  However, 
the  characteristic  frequencies  are  not  the  same  and,  consequently, 
neither  are  the  physical  frequencies  for  which  the  adm.ittances  are  equal. 

The  characteristic  frequency  of  noncircular  lines  will  be  higher 
than  that  of  circular  lines  of  the  samie  cross-sectional  area.  For  any 
excitation  frequency  ui,  this  implies  that  the  characteristic  admittance 
(impedarice)  of  noncircular  lines  will  be  lower  (higher)  than  for 
circular  lines  of  the  same  area,  and  the  designer  is  unfortunately 
faced  with  the  resultant  mismatch.  Nevertheless,  the  characteristic 
fr'equency  (or  radi-os ) of  the  noncircular  line  together  with  the  circular 
line  theory  may  be  used  in  a simple  procedure  to  predict  the  performance 
of  the  noncircular  line  at  any  desired  frequency. 
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The  ratios  uc-p/ij^cc  “Cp/'^Cc  plotted  in  Figures  25  and  27 
for  rectangular  and  annular  lines,  respectively,  assuming  ecual  flow 
areas.  Similar  curves  could  be  constructed  for  any  noncircular  line. 

The  characteristic  frequency'  of  circular  lines  is  easily  determined 
(ijuc^  = St^v/Ac)  and  from  curves  sim.ilar  to  Figures  26  and  27  the  char'a- 
cteristic  frequency  (wc)  of  the  particular  noncircular  line  is 
determined.  For  the  physical  frequency  of  interest,  the  frequency 
ratio  is  computed  and  any  desired  performance  parameter  such  as 

a/..  Figure  IS,  or  |Yo/''^soI>  Figure  24,  is  easily  obtained  from  the 
circular  line  performance  curves.  It  is  emphasized  that  laminar  mean 
flow  is  assumed;  combiined  turbulent  mean  flow  and  laminar  oscillatory 
flow  is  discussed  in  Section  V.  Furthermore,  no  provision  is  made  in 
any  of  the  models  to  adequately  account  for  abrupt  or  discontinuous 
changes  in  transmission  line  flow  area,  or  the  generation  and  propagation 
of  vortices  or  swirl  areas  commonly  found  at  these  discontinuities. 
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Figure  26.  Ratio  of  Characteristic  Frequencies  of  Rectangular 
to  Circular  Lines 
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IV.  The  Trar.sient  Response  of  Koncircular  Fluid 
Transr.issior.  Lines 

The  transient  response  of  fluid  transrriission  lines  of  circular 
cross  section  was  oltained  by  Brown  [2],  who  derived  the  propagation 
operator  and  characteristic  irr.pedarice  of  aii  infinitely  long  uniforrr. 
rigid  line  in  the  Laplace  domain  for  both  compressible  and  incompress- 
ib'le  fluids.  Brown's  solution  is  identical  to  that  of  Iberall  [1]  and 
Kichcls  [3]  frequency  domiain  models  for  the  special  case  of  sinusoidal 
excitation  of  the  circular  line.  Subsequently,  Karam  [19,2C]  developed 
a simplified  solution  for  the  transient  response  of  circular  lines 
after  transforming  approximate  high  and  low  frequency  domain  models 
into  the  time  domain  and  then  observing  a basic  similarity  between  the 
transformed  solutions.  Karam's  solution  is  more  easily  applied  than 
that  of  Brown  [2],  and  may  be  used  to  obtain  the  impulse  and  step 
responses  of  semi-infinite  fluid-filled  circular  lines  with  good 
accuracy  and  without  the  aid  of  a computer. 

In  Section  III,  a method  of  expressing  the  frequency  dependence  of 
noncircular  lines  entirely  in  terms  of  equivalent  circular  lines  was 
presented,  and  results  were  obtained  for  rectangular  and  annular  lines 
over  a wide  range  of  aspect  and  radius  ratio.  The  primary  purpose  of 
this  section  is  to  present  a method  of  transforming  the  frequency 
response  of  noncircular  lines  into  the  time  domain.  The  analytical 
approach  essentially  combines  the  approaches  of  [2,3,20]  and  that  of 
the  previous  section,  and  results  in  a simple  but  complete  solution  of 
the  step  response  of  noncircuuar-  lines. 


Frecuencv  Dciriain  Eauations 


The  exact  analytical  expressions  for  the  parameters  Z,  Y,  T,  and 
Yq  ^hich  characterize  the  perfonr.ance  of  a fluid  transrr.issicn  line  are 
relatively  curherscme  tc  use  in  the  frequency  or  Laplace  domain. 
Furthermore,  it  has  been  observed  that  the  timie  transformis  associated 
with  these  parameters  have  not  yielded  closed  form  analytical  expressions 
ever  any  of  the  frequency  or  Laplace  domain,  even  for  the  simiple  case 
cf  circular  transm.ission  lines.  The  mathemiatical  difficulties  may  be 
traced  directly  to  the  underlying  frequency  or  Laplaoe  domain  solutions, 
which  contain  ratios  cf  Eessel  functions.  These  ratios  are  not  thus 
far  expressible  in  a font  amenable  to  transformiation  into  the  time 
dom.ain.  Furtner,  according  to  Sneddon  [21],  it  is  unlikely  that  the 
Bessel  solutions  will  ever  be  totally  transformed.  Instead,  high  and 
low  frequency  approximations  are  made,  and  these  approximations  are  then 
transfonried  into  the  time  domain.  In  the  mid- frequency  range,  very 
complicated  numerical  techniques  have  been  used  to  accomplish  the 
transfonijation  [22].  Alternately,  the  transient  response  in  the  m.id- 
frequency  range  may  be  approximated  in  the  time  domain  itself,  as  was 
done  by  Karami  [20]. 

The  general  approach  taken  herein  will  be  to  partition  the 
frequency  domiain  into  three  regions,  to  analyze  the  equivalent 
electrical  circuit  of  a general  transmission  line  in  each  of  these 
regions,  arid  then  to  express  the  results  in  the  time  domain.  In  this 
approach,  the  text  by  Weber  [11]  is  especially  useful. 

The  three  frequency  regions  are  defined  with  reference  to  the 
characteristic  frequency,  wc  = Ry/Lv,  ad  follows:  low  frequency 
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(u)  <<  Uq) . Tic-frequency  (w  = w^),  and  high  frequency  (w  >>  Wq).  The 
characteristic  frequency  ratio  = uLv/?'v  be  thought  of  as  a 

frequency  dependent  (AC)  Reynolds  number.  At  low  frequency  ratios 
viscous  flow  effects  are  prevalent  while  at  high  frecuencv  ratios 
inertial  effects  prevail.  In  the  mid-frequency  range  viscous  and 
inertial  effects  are  of  the  sarrie  order  of  magnitude.  Similarly,  the 
thermal  characteristic  frequency  u-  = uc/o^  roughly  divides  the  heat 
transfer  effecrs  into  regions  of  isothermal  behavior,  u « wm  and 
adiabatic  behavior,  u >>  L>y.  Since  the  Prandtl  number  (c^)  is  near 
unity  for  air  (and  man}^  other  gases)  either  u)q  or  up  will  suffice  as 
the  overall  behavioral  index,  and  u^;  has  been  chosen  for  this  work. 

For  liquids  (y  = 1),  there  is  essentially  no  dependence  on  the  Prandtl 
number,  since  terms  in  the  energy  equation  containing  the  Prandtl 
number  are  multiplied  by  the  quantity’  (l-y).  [See,  for  example, 
equation  (26).] 

At  low  frequencies,  the  equivalent  circuit.  Figure  2,  is  comprised 
of  essentially  constant  value  elements  of  R,  L,  and  C with  G = 0. 

These  values  are  R = Ry,  L = Lv>  and  C = = yCa*  -or  any  line  of 

constant  axial  cross  section.  Using  these  values  in  the  defining 
equations  of  the  series  impedance  and  shunt  admittance  [equations  (28) 
and  (2S)],  together  with  equations  (32)  and  (34),  the  low  frequency 
characteristic  impedance  and  propagation  operator  are  easily  determined: 


Yt 


(88) 


= N 2l^L  ' = + IBl  = ( 


yR,/Cqu 


T 


(l+i) 


(89) 
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Therefore,  the  lov  frequency  attenuatior,  or^,  and  phase  angle  increase, 

Bl> 


i-  . /fiT?  . ^ 


WO)- 


(90) 

(91) 


At  high  frequencies,  the  induc‘"ance  and  capacitance  approach  their 
adiabatic  values.  However,  the  resistance  and  conductance  are  frequency 
dependent  and  becorie  large  with  frequency.  Karam  and  Franke  [9]  have 
developed  approxiroations  for  the  high  frequency  circuit  components, 
which  are  rewritten  in  terms  of  as  follows: 


% " "I  L'a  A ' 


Lh  = La  = p‘'7a 


G--:  = Cs  /\|  XluWc  ' 


- 1~- 
'n  " 2o 


Ch  = Cg  = A/yp'" 


(92) 

(93) 

(94) 

(95) 


In  the  high  frequency  region  R <<  wL  and  G <<  wC.  Using  these 
inequalities,  the  high  frequency  attenuation  per  unit  lengtn,  aj^,  and 
phase  angle  increase  per  unit  length,  6^5  approximately  [11]: 


LhSh  + Rh^k 


2‘ uhCk 


(96) 


(97) 
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Substituting  equations  (92) 


through  (95)  into  equations  (96)  and  (97), 


= 4 


v'P  o 

a 

V = 4ca 

1 

c 

r] 


(9S) 


Sh  = ^ 


(99) 


It  is  important  to  note  that  the  foregoing  development  is  applicable  to 
any  fluid  transmission  line  of  arbitrarj,’  but  constant  axial  cross 
section. 

From  equations  (90)  and  (96),  it  is  seen  that  both  and  ajj  ar^e 
of  the  form  a = Kvar,  where  K is  constant,  and  such  a line  has  a rather 
simple  time  domain  solution. 


Time  Domain  Equations 

The  characteristic  time  tj;  is  defined  as  the  inverse  of  the 
characteristic  frequency  (tc  = l/wj,),  and  the  time  domain  is  divided 
into  three  regions  corresponding  to  the  previously  defined  frequency 
regions.  Thus,  the  short  time  region  (t  « t^)  corresponds  to  the 
high  frequency  region  (w  >>  wq)  and  the  long  time  region  (t  >>  tc) 
corresponds  tc  the  low  frequency  region  (u  « Uc). 

Metzger  and  Vabre  [23]  have  shown  that  when  the  high  frequency 
attenuation  is  of  the  form  a.yf  = k/T,  the  step  response  tc  a step 
input  U(t)  at  £ inches  down  a semi- infinite  line  is 


^f(t-T)^ 


1/21 


U(t-T) 


(100) 
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wr.ere  ‘J(:)  * u(:)  is  a ur.is  st*;  ir.pus,  T * i/c  j is  the 

adiaiatic  delay  ti!te,  ar.d  f is  the  sigr.al  frequency.  Tzr  long  tires, 
.^ara.t.  Crc]  has  shown  that  the  unit  step  reipcr.se  say  he  written  as 


Since  fcr  long  times  t » t-T,  evidently  the  'unit  step  responses  for 
long  and  short  tis.es  are  gcverr.ed  hy  the  sane  relation.  It  seems 
reascr.ahle  tc  assure  that  the  step  response  will  retain  the  same  fort 
in  the  si d- frequency  ra.nge.  Thus,  in  general 


(102) 


I 

I 


t 


It  is  noted  that  this  assumption  avoids  the  transfcrmation  of  the  mid- 
frequency  region  into  the  time  domain.  Browr.  and  !:elscr.  [22]  achieved 
some  success  in  the  mid-frequency  range  for  liquid-filled  circular 
lines,  but  the  method  was  numerical  and  was  reportedly  quite  complex  and 
sensitive.  The  assuitiption  that  the  attenuation  is  of  the  form  a = k/T 
in  the  mid- frequency  range  is  a convenience  which  allows  the  entire 
time  domain  solution  to  be  expressed  as  a single  complementary;  error 
function.  The  "constant"  K is  actually  frequency- dependent  in  the  mid- 
range. Karair.  [20]  used  the  time  transform  of  the  known  low  and  high 
frequency  asyrrjptotes  obtained  from  the  equivalent  circuit  approxin.ations 
and  connected  these  with  a logarithmic  straight  line  of  the  form 
Iogio(a/OH)  = A + B logio  where  t'  = tv/ro^.  This  straight  line  patch 
is  therefore  made  in  the  time  domain  and  any  correspondence  with  the 
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-o-air.  is  Icsr  in  th«  nid- frequency  range.  Karar.  apprcxiaated 
the  transient  tehavicr  of  a/au  far  circular  lines  as  [20]: 

1 for  t'<  0.0125 

KC-’T*  for  0.C125  < t < 1.25 

— — r for  t'>  1.25 

1 t ^ 

vhere 

f -3610 


1 t-  ^ 


and  K = (0.0125) 


-c 


A qualitative  plot  of  a/ojj  versus  ^ is  given  in  Figure  26,  along  with  a 
representation  of  Nichols'  transfcr!r.ed  solution.  It  is  er-phasized  that 
the  dotted  line  cf  Figure  26  representing  Nichols'  solution  is  only 
qualitative;  if  the  exact  transfcrrt  of  Nichols'  solution  were  known, 
there  would  be  no  need  for  the  straight  line  approxir.aticn.  Karan, 
defines  tq  = and  expresses  the  unit  step  response  as 


(lOU) 


This  solution  gives  one  of  the  simplest  representations  of  the  step 
response  of  circular  lines  arid  shows  very  good  agreement  with  the  much  mere 
complicated  emalytical  and  n-jr:erical  solution  of  Brown  and  Nelson  [22] 
for  liquid  lines.  Kara.T.'s  solution  is  also  in  good  agreement  with  his 
own  experiii.ental  data  arid  that  of  Kantola  [2*4]  for  air-filled  lines. 

The  pressure  step  response  of  liquid-filled  lines  obtained  using  the 
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Figure  28.  Logarithmic  Approximation  of  Attenuation  Ratio  for  Circular  liiii**s 


results  cf  Ercwr.  e:.d  .‘.'elstr.  1 22]  and  of  K.arar.  [2j]  are  plotted  in 
Figure  25  for  selected  values  of  Karas' s results  are  applicable 

tc  liquids  when  y is  unity  and  the  fluid  bulV  soduluE  Ki  is  used  tc 
cor.pute  the  celerity.  In  [22, 2^*],  the  ncndir.ensicnal  tir.e  is  also 
defined  as  t'=  tv/r^*.  It  is  seen  fjxr  Fif-re  29  that  the  predicted 
responses  are  in  good  agreerrent  over  a wide  ranige  of  to- 

Pressure  (Flow)  r.estonse  of  Noncircular  lines  tc  a 
Step  In; ut  cf  Pressure  ( Flow ) 

Whereas  Karar.  and  others  used  the  quantify  v/rj,^  tc  ncndiaension- 
alize  tine,  the  characteristic  tine,  tj.  = will  be  used  herein. 

Ihius,  t' and  tV  are  redefined  and  a characteristic  tine  equal  to  unity 
is  introduced  as  fellows: 

■m  ■ 

T = tu~Z  = t/tc  ■ 


TC  = U^Z  = T/tj 


6Tq 


(1C5) 


1 = Wgtc  = *c/^C 


In  terns  cf  the  redefined  nondimer.sior.al  tines,  equation  (10**)  becomes 


Pj.  = erfc 


• ^1  \k] 


U(t-To) 


(loe) 


and  the  expression  for  the  long  and  short  tine  step  responses  become 


2 V T-Tr 


1/2 


U(t-to) 


(107) 
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•r:-;  •*-:  SCT-Tc) 


(i:6) 


The  Icgarithnic  straight  lir.e  ap p rex it;a tier.  r.a-/  he  rewritter.  as  fellows; 


I 

for 

T < 

0.1 

Kt'^- 

for 

0.1 

< t < 10 

.-ey 

y-1 

for 

T > 

10 

(IC'5) 


where 


E = -r  logio 


y-1 


and  K = (0.1) 


Equations  (106)  and  (109)  approximately  describe  the  step  response 
of  lines  of  constant  axial  cross  section,  either  circular  or  noncircular. 
The  key  point  is  the  definition  and  use  of  since  the  ratio  w/u;^ 
(hence,  t = t/tj,)  characterizes  the  response  of  the  selected  transmission 
line.  For  long  times,  say  one  decade  or  more  above  t^.  (t  i 10),  the 
line  response  is  characterized  by  predominantly  viscous,  isothermal 
behavior.  At  one  decade  below  tj.  (t  S 0.1)  the  line  performance  is 
nearly  inertial  and  adiabatic.  Furthermore,  the  step  response  of  all 
lines  of  the  sairie  nondimens ional  delay  time  tq  = ui^T  may  be  represented 
by  a single  curve  when  plotted  against  t = w^t . Equation  (106)  versus 
t is  plotted  in  Figure  30  for  several  values  of  tq.  The  fluid  has  been 
assumed  tc  be  air  at  80  F (y  = 1.4017).  Figure  30  represents  either 
the  pressure  response  downstream  to  an  upstreanj  pressure  step  input  or 


81 


I 


o 


Figure  30.  Pressure  (Flow)  Respons#*  to  a St«»p  Input  of  Pressuro  (Flow)  for  Air-Flll«*cl 
Liner,  of  Arbitrary  Cross  Section 


he  flow  resrcr.se  .icwr.srreax  ro  er. 


:..e 


T.'.e 


-pstres.r.  f.cw  step  ir.put, 
equivalent  transfer  furctior.  ir.  the  Laplace  cor.air.  is 
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Fress’ure  (.^low)  Kesttr.se  ^ Nontircular  Lines  a 
Step  Input  c_f  Flow  Pressure) 

Twc  other  transfer  functions  are  of  interest  in  deterrining  the 
response  of  ser.i-ir.f inite  lines.  These  are:  (1)  the  pressure  response 
to  a step  input  of  flow,  and  (2)  the  flow  response  to  a step  input  of 
pressure.  The  respective  Laplace  transfer  functions  are: 


?r<s)  _ 1 , , ^'ir(s) 


(111) 


Qr(s) 

Vl7 


iYo(s)  e 


tr(s) 


(112) 


and  it  is  noted  that  a convolution  integral  nay  be  used  to  solve  for  the 
time  domain  response.  Karar.  [2C]  has  used  the  following  equations  to 
approximate  these  responses; 

p ft  p p ft 

7^  (t)  = Lo(A)  (T-A)dA  ^ (t)  Zo(l)dA  (113) 

"s  Jq  ''s  Jq 

r ft  p ■p  r^ 

(t)  = I YoCA)  ^ (T-A)dl  ? ^ (t)  I Yo(X)dX  (114) 

s Jq  ;o 
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(>.)dX  = L 


-[¥] 


Yo..(X)d>.  = L' 


For  t'  s ■z\/Tq-  and  t = u-t  equations  (117)  and  (115)  Lecone 


(116) 


(115) 


(120 


I 2c.(>)dX.Z, 


(X)dA  = Y. 


■ ^ -K’  T 


(121) 


(122) 


It  is  convenient  and  acre  realistic  to  construct  a single  continuous 
f'jTiCtion  to  describe  the  behavior  of  the  integrals  of  Zq  and  Y^. 

Karair  [20]  has  foiu^d  that  equations  (11£)  and  (121)  nay  be  combined  as 
follows  to  give  a closer  approximation  of  the  overall  behavior  cf  the 
integral  of  Z^. 


(123) 


Similarly  a continuous  function  approximating  the  integral  of  Y^  is 
obtained  in  terms  of  t ' and  t as 
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(12-) 


Equations  (125)  ar.c  (12- ) versus  t'  are  plo-.t«d  ir.  figure  31,  al&r.g  with 
the  asyrtctotic  lin.its  giver,  t'v  equaticr.s  (116),  (120),  (121),  ar.i  (122), 
for  air-filled  circular  lir.es  with  v = l.i*tl7.  It  is  seer,  that  t.he 
twc  ccr.tinuous  functicns  of  equations  (123)  and  (12'«)  are  properly 
bounded  by  the  respective  asycptctes.  Using  equations  (113),  (11*«), 
(123),  and  (lO**)  the  desired  step  responses  are  obtained  ir.  dittension- 
less  fcrrr.  as 


(125) 


-1/: 


(125) 


Equations  (125)  and  (126)  are  plotted  in  Figures  32  and  33,  respectively, 
along  with  the  results  of  Brown  and  Nelson  [22]  for  liquid-filled  lines 
(•V  = 1)  and  t'  = tv/r^^.  It  is  seen  that  the  two  approaches  are  in 
good  agreement  over  a practical  range  of  t ' . The  above  step  responses 
are  adapted  to  non-circular  air-filled  lines  as  was  done  for  Fr'^s* 
the  results  are  shown  in  Figures  3“  and  35.  Ir.  these  latter  figures, 
the  dimensionless  times  are  based  or.  the  characteristic  frequency, 
i.e. , t = w^t. 
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33.  Comparison  of  Flow  Response  to  a Step  Input  of  Pressure 
Circular  Lines 


Figure  34.  Pressure  Response  to  a Step  Input  of  Flow  for  Air-Filled  I,in‘-s  of  Arl'itr'ary 
Cross  Section 


cross  section  cf  higher  In  other  words,  lines  of  the  same  tq  and 

working  fluid  properties  but  different  Uq  necessarily  imply  different 
axial  locations  (lengths)  along  the  semi-infinite  line.  It  follows 
that  the  distance  along  the  line  will  be  shorter  for  the  line  cf  higher 


The  line  of  higher  will  also  attenuate  high  frequencies  at  a 
mere  rapid  rate,  and  will  approach  the  behavior  of  the  diffusive  RC 
line  more  rapidly.  For  signal  transmission,  lines  of  higher  should 
be  avoided,  as  these  are  accompanied  by  larger  attenuation  per  line 
wavelength  at  any  physical  frequency. 

As  tq  becomes  large,  say  tq  2 1,  the  line  becomes  essentially  a 
diffusive  RC  line,  with  R = R^  and  C - C^.  For  a circular  line  of 
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ccrrespir.iir.g  length  is  122  ft  and  7 is  0.10'’  sec. 


Fcr  the  sa~e  line  cjcss  section  and  t;  . rut  filled  >.-it".  water  at  E'C®:, 
the  corrcstcncing  length  is  ejS'c  feet  and  T is  l.°4i  seconds.  It  is 
unlikely  that  signal  transmission  would  be  cf  interest  ever  such  long 


fluid  lines.  For  such  lines,  treatment  as  a lumped  ?.C  line  would  be 
sufficient  fer  most  purposes. 


Sumcc.ary  of  Transient  responses 

In  this  section  a method  was  presented  for  determ.ining  the 
transient  response  of  semi-infinite  fluid-filled  line: of  arbitrary’  but 
constant  axial  cross  section.  The  results  were  presented  in  terms  cf  a 
characteristic  frequency  which  depends  only  on  the  lam.inar  steady  flow 
resistance  and  conductance  cf  the  arbitrary  line.  Equivalent  circuit 
paramjeters  (R,  G,  L and  C)  were  developed  and  used  to  obtain  the  time 
transfonr.s  cf  the  various  combinations  of  pressure  and  flow  responses 
to  step  inputs  cf  pressure  and  flow.  Results  were  compared  with  these 
cf  ethers  for  both  compressible  and  incompressible  fluids.  For  non- 
circular  lines  the  results  are  shown  for  air  at  80®F.  It  is  noted 
that  the  results  are  restricted  to  la.t.inar  flow.  Similar  results  for 
the  various  step  inputs  propagating  through  fully  developed  turbulent 
miean  flow  are  presented  in  Section  VI. 
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La-.inar  Cscillatory  Zort  ir.ed  vith  Steady 

Zeveloted  Mean  “lov 


Sac>:ground 

Several  irivestigatior.s  have  been  corrieted  concerning  the  trans- 
mission or  small  oscillato3rv  sigTials  through  fully  developed  mean  flow 
A majoritv’  of  these  have  been  directed  toward  the  superposition  of 
laminar  oscillatory  signals  on  fully  developed  laminar  mean  flow  [2,3, 
Recently,  investigations  concerning  the  behavior  of  oscillator^.’  signal 
combdned  with  fully  developed  furbulent  flaw  have  been  comoleted  [li, 
25,27],  Most  of  the  recent  efforts  have  been  restricted  to  inoomoress 


adiabatic  flow  in  circular  lines. 

Kichcls  [3]  obtained  the  small-signal  frecuencj'  response  of  fluid 
filled  circular  lines  including  frequency  dependent  dispersion  of  the 
propagation  speed  and  distortion  of  the  oscillatory  signal  wave  shape. 
In  addition,  heat  transfer  effects  were  included.  In  essence,  .N’ichols 
separated  the  steady  and  oscillatory  fluid  parameters  (as  was  done  in 
Appendix  C)  and  subsequently  solved  the  oscillatory  flow  problem; 
nevertheless,  his  solution  is  applicable  to  flows  combining  steady 
developed  laiijinar  mean  flow  (Kagen-Fciseuille  flow)  and  laminar 
oscillatory  flow.  HcLmboe  and  Rouleau  [26]  have  shown  that  if  the  ini 
conditions  are  assumed  to  be  .Hagen-Poiseuille  flow  then  the  governing 
differential  equations  may  be  placed  in  a form  which  is  identical  tc 
Kichols'  equations.  The  results  of  Brown  [2],  expressed  in  both  the 
Laplace  and  time  domains,  are  identical  tc  those  of  Kichols  for  the 
special  case  of  sinusoidal  excitation  signals. 
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t.'.e  actTiOrs'  ricdeis.  F.ecer.tlv  Fur.k  ar.d  Wocc  [2'j  develcoec  a sir.gle 
bo-ndaiy  layer  model  to  accour.t  for  the  trar.sier.t  viscous  effects  ard 
applied  this  node,  to  snail  ar.plituce  sir/jsc idally  disturbed  turbulent 
flew  in  circular  lines.  The  authors  developed  a transfer  function 
relating  the  local  boundary.'  layer  velocity  gradient  to  the  core  velocity 
of  the  flow.  Tne  model  is  comparatively  much  simipler  than  that  of  Browr. , 
e_t  £l . [15];  nevertheless  it  gives  values  cf  frequency  dependent  signal 
attenuation  which  are  in  excellent  agreement  with  the  theory  arid  the 
experiiTiental  data  of  Reference  [15].  It  is  noted  that  the  utility  cf  the 
riodel  given  in  Reference  [25]  is  limiited  to  the  higr.  frequency  regime  and 
constant  phase  velooity;  .'"unk  and  Wood  suggest  the  use  of  a quasi-steady 
model  (i.e.  , a constant  IRC  mjodel)  at  low  frequencies.  It  will  be 
subsequently  shown  that  the  model  of  Funk  and  Wood  is  closely  related  to 
the  high  frequency  approximations  developed  by  Nichols  [3]  and  Karam.  [2C] 
for  laminar  oscillatory  flow  cciTibined  with  laminar  steady  flow,  and  to  the 
concept  of  a characteristic  frequency  and  radius  previously  developed 
herein. 

Trikha  [29]  has  suggested  that  the  frequency  dependent  part  of  the 
friction  in  turbulent  flow  may  be  approximated  as  the  frequency  dependent 
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•art  or  the  mctror 


larir-ar  flew  for  a lir.ited  rar.ge  of  r.eynclcs 


r.a.-Tler  (in  this  r.atter,  see  also  [7]). 

in  this  section  a siir.tlifiec  Ttodel  of  oorrilinec:  steadv  turbulent  and 

tne  cor.tcnents  ?,,  G,  1,  and  C of  the  eo_uivalent  electrical  circuit, 
rigure  C,  and  also  or.  the  concert  of  a characteristic  frequency,  h 
nethod  cf  detenr.ining  the  frequency  response  of  lines  of  either  circular 
or  ncr.circular  cross  section  carrying  fully  develctec  turbulent  ji.ean  flow 
is  developed.  The  results  are  presented  fer  a dic.ensionless  attenuaticn 
and  frequency.  For  the  circular  line  rodel,  results  are  compared  with 
these  given  in  [15,25]. 

In  Section  VI,  Karam's  approximate  method  [19]  is  adapted  to  lines 
of  both  circular  and  ncncircular  cross  section  carrying  turbulent  mean 
flow  and  the  pressure  and  flow  transient  responses  cf  the  lines  to  step 
inputs  are  estimated.  Before  proceeding  with  the  broader  purposes  of 
this  section,  it  is  instructive  to  review  the  physical  aspects  of 
comi'ined  oscillatory  flow  and  steady,  developed  mean  flow  cf  compressible 
or  incompressible  fluids  in  a semii- infinite  circular  line. 

Physical  Considerations  of  Comb,ined  Oscillatory  and 
Steady  Developed  Mean  Flow 

A large  body  of  literature  exists  concerning  steady,  fully  developed, 
laminar  or  turbulent  flow  in  circular  lines,  from  which  one  may  select 
flow  relationships  ranging  in  complexity  from  relatively  straightforward 
to  very  complex.  For  turbulent  mean  flow,  elegant  theoretical  flow 
descriptions  involving  statistical  approaches  and  correlation  techniques 
have  been  proposed.  However,  these  descriptions  have  failed  to  yield 
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sclutior.s  except  ir.  the  sir.plest  cases.  The  ap 
tc  develtp  the  Ciecessary  steady  flov  pararr.eters 
pher.rreT-.rlogical  descriptions  of  steady  turbule 
spprcach,  the  texts  cf  Schlichtir.g  [5C]  and  Hin 
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For  stall  signal  excitation  of  developed  turbuien 
flow  velocity  profile  will  not  be  greatly  altered  by  t 
oscillatory  signal  provided  the  steady  flow  Feynclds  n 
cier.tly  large,  and  the  cscillatcry  signal  an-.plifcde  is 


,t  flow,  the  steady 
he  superitpesed 


urber  is  suffi- 
rfc-»dtivej.v  sr.a.,,*. 


Steady  and  Oscillatory  Boundary  Layers 

Fcr  steady  developed  latiinar  flow  in  a circular  lire  or  duct  the 
boundary  layer  is  presumed  to  extend  from  the  wall  tc  the  centerline  cf 
the  circular  cross  section.  Thus  the  steady  flow  laminar  boundary  layer 
thickness  denoted  is  always  equal  to  the  radius  cf  the  circular 

line  = Tq),  and  is  not  dependent  on  the  laminar  flow  .Reynolds 

numiber.  The  classical  .Hagen-Foiseuille  velocity  profile  is  obtained 
for  all  laminar  flow  in  circular  lines.  However,  for  steady  developed 
turbulent  flow  the  associated  velocity  profile  and  boundary*  layer  thick- 
ness will  depend  on  the  Reynolds  number.  As  the  t'orbulent  flow  Reynclds 
nurrb>er  increases,  the  associated  velocity  profile  changes  shape  and 
approaches  a plug-flow  profile.  Simultaneously  the  associated  turbulent 
bour.daiy  layer  thickness  decreases,  indicating  large  velocity 
gradients  at  or  near  the  wall  and  large  resistive  losses. 

For  laminar  oscillatory  flow,  the  dynamic  boundary  layer  thickness 
will  vary  with  frequency.  Nichols  [3]  proposed  a boundary  layer  thick- 
ness, 6^  “ t^2v/u),  fcr  circular  lines;  and  it  is  seen  that  the  oscillatory 
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bc\;rdar\’  layer  thick-ess  decreases  with  ircreasir.g  .'rec/jency . Si-ilar 
effects  ccc'cr  ir.  electrical  conductors  where  is  called  the  skir.  deoth. 
At  lew  frecyuer.cies  the  cscillatcrn'  velocity  prefile  attrcaches  the 
classical  Hagen-Pciseuille  cistribution. 

■.'.’hen  la-T.inar  oscillatory  flow  is  corb-ined  with  steady  rear,  flow, 
it  is  assured  that  an  interaction  between  the  steady  flow  boundary  layer 
anc  the  dynaric  (oscillator;.’)  boundarv  laver  will  occur.  If  the  steady 


flew 

boundar.’  layer  thi 

cVd^ss 

f s 1 

ess 

than  that  of  the  oscillatory  flow 

(i . e. 

‘5  C'vc  ^ '-'i  OT  Cv  — 

< fi,) 

then 

the 

oscillatcry  flc’w  will  be  derenier.t 

uecTi 

the  steady  flow  pr 

'cfile. 

On 

the 

other  hand,  when  (or 

< 

iv-)  the  ’oehavior 

of  the 

osci 

Ilf 

torv  flow  tends  to  become  indetendent 

cf  the  rear.  flew.  This  observation  is  substantiated  both  analytically 
and  experirentally  by  the  work  of  Brown,  et  al.  [15],  "ranke,  Karan,  and 
Lirreumer  [1^*],  and  others  [9,25].  Stated  in  other  terrs,  the  behav’ior 
of  fne  oscillatory  signal  in  the  low  frequency  region  depends  on  the 
rean  flow  -whereas  in  the  high  frequency  region  the  behavior  of  the 
oscillator)'  signal  becomes  essentially  independent  of  the  mean  flow  and 
dependent  on  the  oscillatory  signal  frequency  (x).  For  the  experimental 
data  available,  the  high  freq’uency  behavior  appears  to  approach  laminar 
oscillatorv’  flow  behavior  [15,27].  It  is  f'urther  assur.ed  that  for  an’y 
steady  flow  there  exists  a frequency  at  which  the  dynaniic  and  steady  flow 
boundary  layer  thicknesses  are  comparable  (6^  ^ 6v).  In  this  region,  a 
transition  from  essentially  steady  flow  dependence  to  frequency 
dependence  takes  place,  or  vice  versa.  Similar  qualitative  behavior 
would  be  expected  concerning  the  temperature  profiles,  particularly  for 
fluids  with  a Frandtl  number  near  unity.  For  circular  lines  Nichols  [3] 
defines  a thermal  nonadiabatlc  skin  depth  as  67  = /2V7/W  which  is  seen 


98 


tc  be  frequency  ceyenient.  A.t  lev;  frequencies  the  heat  transfer  behavior 
cf  the  oscillatory  flow  will  depend  essentially  or.  the  rear;  flow, 
whereas  at  high  frequencies  the  heat  transfer  will  depend  essentially 
upon  the  oscillator-/  frec/uencr.- . y-uch  cf  the  discussion  cf  the  frequency 
dor.ain  equations  in  iecticr.  IV  applies  qualitatively  to  the  al.cve 
discussion.  It  is  net  intended  to  extend  the  discussion  of  the  various 
boundar/v-  layers  any  f-urther,  since  ultiirately  a Siodel  for  circular  and 
noncircular  lines  with  rear,  turbulent  flow  is  sought.  For  certain  non- 
circular  lines,  ncre  than  one  or  fv  -s  involved  and  it  is  rore 
convenient  to  eliminate  the  dependence  of  the  flow  cr.  the  georetric 
cross  section  of  the  line.  By  integrating  the  various  parameters  over 
the  appropriate  cross  section,  geor.etrq.-  dependence  may  be  essentially 
removed,  as  was  previously  done  in  developing  the  characteristic  frequency 
and  components  cf  the  equi->^alerjt  electrical  circuit  for  la.minar 
oscillatory  flow. 


A Mean  Flow  Characteristic  Frequency  for  Circular  Lines 

The  concept  of  a characteristic  frequency  -which  divides  the  flow 
into  regions  of  viscous,  isothermal  behavior,  and  inertial,  adiabatic 
behavior  has  been  discussed  in  Section  III  for  circular  and  noncircuiar 
lines.  It  is  desired  to  extend  this  concept  to  lines  of  arbitrain- 
cross  section,  carry-ing  combined  oscillatory-  and  steady  turbulent  mean 
flow.  Returning  to  the  equivalent  electrical  circuit,  Fig-ure  2,  and 
keeping  in  mind  that  the  underlying  steady  flow  determines  the  low 
frequency  behavior  of  the  oscillatory  signal,  it  is  necessary-  to  determiine 
R,  L,  G,  and  C for  the  underlying  steady  flow.  The  low  frequency- 
resistance,  Rv,  is  dependent  on  the  mean  flow. 
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and  similarly  for  tr.e  inertar.ce, 


Lv  = Ki 


1 

i^)‘  dA 

A 

Vu/ 

J 

A 

illz) 


It  the  underlying  it, ear.  flow  is  fully  developed  laitinar,  it  was  treviously 
shown  that  for  circular  lines 


Kvj  = 


6ru 


u 

= — L 


where  suhscript  i denotes  la.t.inar  flow  conditions.  These  parair.eters  do 
not  depend  or.  Reynclds  nurjrer.  For  turbulent  tiear.  flew,  the  shape  of  the 
velocity  profile  is  a function  of  the  friction  factor  and  Reynolds 
nuirier.  The  defining  relation  for  the  steady  flow  friction  factor,  f, 
for  lar.inar  or  turbulent  njean  flow  in  circular  or  ncncircular  lines,  is 


-3t  . f 1 -2 
32  ■ D..  T 

fi 


(129) 


rroTT.  equations  (127)  and  (129),  the  lair.inar  or  turbulent  steady  flow 
resistance  of  any  line  of  arbitrary  but  constant  cross  section  is 


(fRe)u 


AIv  2AD.^^ 


(130) 


The  simplest  fonr,  of  tne  relation  between  f and  Re  for  developed  turbulent 
flow  in  both  circular  and  noncircular  lines  is  the  Elasius  relation. 


f = a 3164/Re 


1/4 


(131) 
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wr.^r-^  w.'.e  ?eyr.clc3  nurier,  re  is  defined  in  terT.s  cf  tr.e  r.vdreulic 


ciar^eter : 


(132) 


Zquetisr.  (131)  gives  good  accuracy  up  tc  ?.e  < E>1C-  [33],  and  is  used  tc 
ccrcpuce  f r.erein.  Using  Equations  (13C)  and  (131)  the  turbulent  viscous 
resistance  cf  lines  cf  arbitrary  cross  section,  becor/es 


C.  3ie‘*u  ^Z/u 
r-  (Ke) 


(133) 


fror  Hinze  [31],  the  steady  turbulent  velocity  profile  in  circular  lines 
c.ay  be  expressed  as  a siniple  power  law. 


^ 

'^ir.ax  V^o/ 


(13*-) 


where  n is  a function  of  F.e , y = r^  - r,  and  u^^.^x  is  the  r.axir.ui:. 
(centerline)  velocity.  Further,  the  work  cf  Nunner  as  presented  by 
Kinze  [31]  shows  that  the  exponent  n r.ay  be  expressec  as  a function  of 
f only. 


- = vT' 

n 


(125) 


Equation  (135)  is  shown  to  be  in  good  agreement  with  the  experimental 
data  of  Nunner,  Nikuradse,  and  Laufer,  for  both  smocth  and  rough  circular 
lines  [31].  From  Schlichting  [30], 


2n^ 

)(2n+l) 


(136) 
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and  K;,.  is  seer,  fc  be  a sir.ple  expressicr. . In  terr.s  cf  tne  friction 
factor 


^Lr  = 


^ (>'f  t DC*-?  t 2)^ 

^(2/f  + 1) 


(139) 


The  value  of  Kt^-^  will  be  near  unity  since  the  steady  turbulent  velocity 
profiles  are  siir.ilar  to  inertial  or  plug  flew  profiles.  As  ?€  becories 
large,  it  car.  be  seen  that  u/u  and  approach  unity,  while  n - » and 
f -►  0.  It  is  further  noted  that  Ki-j  near  unity  iir.plies  that  bv^  does 
not  differ  much  from  Lg,  and  the  assumption  that  tr.c  frequency  dependent 
inertia  is  in  fact  constant  (L  = La*  or  L = Lv* ) is  a good  approximation, 
particularly  at  high  Reynolds  numbers.  The  appreximatien  Lv»  = la  is 
particularly  useful  in  dealing  with  lines  of  arbitrary  cross  section 
carrying  turbulent  through  flow,  where  expressions  for  u/u  are  either 
very  complex  or  unknown. 

A broad  range  of  values  of  for  circular  lines,  together  with 
closely  related  paramieters , were  obtained  using  equations  (131),  (135) 
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includ*:  ir,  Taile  VI.  It  is  seen  that  the  tahulatec 
with  the  p-retetir.g  -iiscussicr.  cf  . 


.r--.*rr.  w ww^rcN. 


TABLi 

~lcw  Ir.ertante 


Re 

fRe 

^ ■ ( 

kl*  i 

! 

c. 
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0.056265 

56.3 

4.2158 

1 

1 

1.0500  ] 

} 5,0CC 

I 

0.037627 

IBE.l 

C 1 = C ^ 

1.0352 

1C,  COO 

0.031640 

316.4 

5.6219 

1.C302 

1 

50,000 

0.C21159 

1,057.9 

6.6747 

1.0211 

100,000 

0.0177S2 

1,779.2 

7.4569 

i.ciec 

500,000 

0.011895 

5,S4S.S 

9.1675 

1.0125 

X y OOu  y ccc 

1 

0.010005 

10,0C5.C 

9.9973 

1.0106 

1 

i 

Characteristic  Frecuency  and  Radius  - Turbulent  Mean  Flew 

Proceeding  as  in  Section  III,  a set  of  characteristic  frequencies 
and  radii  based  on  turbulent  mean  flow  quantities  are  defined  as 


u 


ct 


L 


8v  ' 
“vt 


(140) 


(141) 


where  subscript  t denotes  turbulent  mean  quantities. 

It  is  more  convenient  to  use  the  characteristic  quantities  defined 
in  equation  (l4l),  since  these  quantities  (i.e.,  adiabatic  inertance) 
are  more  easily  determ.ined,  and  are  consistent  with  the  previous 
discussion  and  other  publications  on  the  subject  [3,25].  Both 
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jrbulerjt  fricticr,  factcr  and  reyncids 


anc  uj-  ara  C6i.€:rider.t  on  tr.s  tui 
r.urber.  In  essence,  the  circular  line  carrying  fully  developed  tiorbulent 
flow  is  replaced  Iv  an  equivalent  '"aller  line  cf  radius  r.^;_  or 
carrving  fuilv  developed  lan.inar  flov;,  of  the  S3i:;e  character istic  frecuenc 
Kichcls'  frequency  donain  solutions  [3D,  applicalle  to  circular 
fluid  lines  carrying  either  developed  lan.iner  mean  flow  or  no  mean  flow, 
are  to  be  used  tc  solve  the  nom.alized  frequency-dependent  perfcrn.ance 
of  bcth  circular  and  noncircular  lines  carrying  fully  developed  turbulent 
flow,  and  are  repeated  herein  for  the  convenience  cf  the  reader.  The 
impedance,  admittance,  propagation  operator  and  characteristic  iitpedance/ 


r = /2V  s a t j6 


(144) 


2o  = 


-t/Y  = 


1/Yo 


(145) 


To  use  the  above  solutions  as  stated,  the  radius  r^,  is  replaced  by  one 
of  the  characteristic  radii,  equations  (140)  or  (141),  and  for 
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turi'^lent  fiov  the  rurther  assuKpttot.  vr  = v is  rtade.  This  ass'Jjt.ttion 
irtlies  a ?rar,dtl  r.ufier  cf  unity  which  is  thought  to  he  Ticre  repre- 
sentative of  turbulent  processes  [30,31].  The  adiabatic  inertar.ee  and 
capacitance  (L^,  C^)  are  cor.puted  fror.  the  actual  line  properties, 
ucuatiens  (1^1)  through  (Iho)  are  the  sar.e  equations  which  were  usee 
tc  clcsely  apprcxir.ate  the  nonr.alioed  frequency  dependence  of  both 
circular  and  noncircular  lines , carrp’ing  either  laminar  or  nc  mean  flow 
and  the  results  were  shown  in  “igures  1£,  19,  iu,  and  3:.  These 
ncrr.aliced  results  will  remain  the  same,  when  plotted  against  the 
characteristic  frequency  ratio, 

As  mentioned  in  Section  III,  for  a given  value  of  the  argument  cf 
the  Bessel  functions,  v''u,/v ' r (or  equivalently,  d/d^  or  the 

normalized  frequency  dependence  will  be  the  same  when  plotted  against 
this  arg'umjer.t.  Since  it  now  is  intended  tc  apply  the  above  modeling 
approach  tc  circular  and  ncncircular  lines  carrying  Turbulent  througn 
flow,  it  is  necessary  to  closely  examine  the  range  of  validity  and 
application  of  this  approach.  The  following  development  will  address 
the  lew,  ir.id,  and  high  frequency  domains  relative  tc  the  turbulent 
viscous  characteristic  frequency,  Wy. , and  Nichols'  solutions. 

Low  Frequency  Domain  - Turbulent  Mean  Flow.  For  u <<  a 
quasi-steady  constant  LRC  model  is  the  simplest  approach,  as  suggested 
by  Brown,  [IB]  and  Funk  and  Wood  [25].  The  equivalent  electrica 

circuit  is  comprised  cf  the  constant  value  elements,  R = Rv^> 

L = L5(=  Lv.^),  G = C,  and  C = Cp  = yC^,  The  lew  frequency  series 
impedance  and  shunt  admittance  are 
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"ror^  equations  fl'-l''  ani  (i^^o),  it  is  seer;  that  as  u attroaches  aero, 
r;_*  at-rrcaches  the  steady  turbulent  DC  resistance,  Ry-,  equation  (133). 
rhe  low  frecuencv  rrotafation  ot-erator  and  characteristic  ir.7)edance  are 


r„  = ‘ ^ (i  - 5 ^)]  ' 


Z-  1 

“t  . <=a 
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"or  very  low  frequencies,  w/wy.  i.  the  above  extressions  ir.ay  be 
further  sirctlif ied: 


_ rr“r — — i 1 rt 1 

Tl.  = i jvRv.-a-  = ^ ^ ^ 


if- 


2,  = \ -Zi  . h 

\ jyu)Ca  A u;  ^ 

^0  = 'I 


1 


2ya,v, 


(l-j) 


(150) 


(151) 


where  c — = 2^0= 


Tu  facilitate  discussion  of  the  constant  LP.C  model  and  comparison  with 
recent  work,  the  dimensionless  frequency,  fi,  and  normalized  attenua- 
tion, oBj.,  for  circular  lines  introduced  by  Brown,  at  [15]  will  be 
used.  These  are 
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Using  the  constant  LRC  nocel,  equation  (l^S),  and  the  aicve  pararr.eters , 
it  can  he  seer,  in  figure  3£  that  a principal  effect  cf  increasing  the 
steady  furiuler.t  fep-nclcs  r.-uther  is  to  increase  the  attenuation  at  all 
frequencies,  up  to  the  laTT.inar  reference  curve,  figure  3£  shows  the 
extension  cf  the  constant  IRC  curves  beycnd  the  lat.inar  limit,  and 
disc  indicates  that  fcr  a given  Heynclds  nuniher  a plateau  is  reached, 
after  which  the  l.RC  model  fails  to  predict  further  increases  in 
attenuation. 

Tc  chtair.  the  limiting  value  of  a,  denoted  av,  for  the  constant 
IRC  mode:,  and  large  values  of  uj,  is  first  rewritten  using 
equations  (146)  and  (147),  and  then  expressed  in  tertTiS  of  its  modulus  and 
phase  angle. 
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Using  the  sHicll  angle  approxiir.ation , sinB  = G, 


(1 


7) 


Siir.iiarly,  fcr  the  plateau  phase  lag,  £;/, 


u r- 

— iy  sin 
Ca 


1 


u 


(156) 


Tne  asymptotes  of  value  ciy^  for  each  Reynolds  nuitber  are  also  clotted 
in  Figure  36,  and  it  is  seen  that  ay  is  an  accurate  representation  of 
the  constant  LRC  attenuation  in  the  mid-frequency  range.  Equating  the 
expressions  for  a given  by  equations  (15C)  and  (157),  the  low  frequency 
corner,  or  break  frequency  ratio,  oi/uiv^  = 1/2,  is  obtained. 

To  convert  to  Brown'”'  dimensionless  frecuency,  used  in  Figure  36, 


The  locus  of  break  frequency  ratios  in  terms  of  57  is 


(159) 


(160) 


( 
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whicr.  cari  be  rearranged  tc  give 


Per  incompressible  flow,  > = 1,  and  equation  (1£2)  becomes 


In  terms  of  Brown’s  dimensionless  frequency,  see  Figure  36, 


and  for  incompressible  flow 


(163) 


(16U) 


(165) 


I 


no 


..Te  cr.tire  set  cr  constant  -.-.C  uiocec  attenuation  curves,  ;igure  zi , 
can  be  reduced  exactly  to  a single  curve,  rigure  5’’,  valid  for  a.l 
F.eynolds  nur.bers  and  line  cress  sections  by  using  the  dinensicnless 
frecuency  ratio  dj/cv,.  and  extressing  the  dimensionless  attenuation  as 


t'v,.‘C3 


dee ) 


vhere  a- _ is  the  real  tart  of  F:,.  in  ecuation  (150). 


It  is  seen  from  eo.uations  (150),  (If”),  and  (152),  and  shown  in 
Figure  37,  that  ay,,  and  all  asynittotic  parameters  are  principally 


functions  of  or  n/uv  • Icuation  (16£),  with  a-*  determined  from, 
equation  (14fc)  and  K = is  plotted  in  Figure  37,  together  with  the 
appropriate  asymptotes  and  corner  frequencies.  A second  curve,  which 
has  been  normalized  such  that  the  plateau  value  of  dimensionless 
attenuation  is  unity  (K  = ),  is  also  shewn  in  Figure  37.  It  can 

be  shown  via  considerable  algebra  that  this  curve  is  also  identical 


to  the  ratio  c/cg  (or  obtained  when  using  t: 


ie  constant 


IRC  model . 


y.id-Frecuenc',’  Domain-Turbulent  ^5ean  Flow 


If  is  used  in  place  of  Tq  in  equations  (1^2)  through  (l*^^), 
bicncls'  circular  line  theory  predicts  the  low  frequency  attenuation 
very  accurately’  at  cr  below  the  low  frequency  break  (w/wy^  = 1/2), 
and  can  be  adjusted  to  predict  the  high  frequency  pararrieters  also. 
However,  in.  the  mid-frequency  range,  w =■  ojy^,  Nichols’  theory  does  not 
predict  the  extent  of  the  attenuation  plateau;  in  this  region  the 
constant  LRC  model  is  used  and  the  plateau  extended  apprcximiately  tc 
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nuation  and  Phase  Velocity  in  CifcuiaT’  or  None  iron  lac  l.ines,  i.ainina 
I Plow,  Constant  PKC  Model 


ur.exrected  anc  drar.atic  variations  in  both  attenuation  and  phase 
velocity  were  encountered  in  the  transition  region  between  the  quasi- 
steady (constant  LRC),  and  high  frequency  (lair-inar  oscillatorpv ) regions 
of  flow  in  water-filled  circular  lines  carrying  turbulent  nean  flow. 

This  san.e  behavior  was  also  observed  by  Vining  [32j  in  a recent 
experisiental  investigation  involving  air-filled  rectangular  lines 
carrying  turbulent  ir.ean  flow.  Vining  and  the  author  conducted  a thorough 
search  for  experimental  anomalies  or  errors  to  no  avail.  The  data  was 
repeatable  from  day  to  day,  after  extensive  expeririental  equipment 
replacement  and  recalibration.  A similar  exhaustive  search  was  made  by 
y.argolis  and  Brown,  and  they  were  eventually  persuaded  that  the  experi- 
mental data  was  correct,  details  of  the  cited  behavior  and  discussion 
of  the  possible  causes  are  best  given  in  [27,32].  However,  the 
behavior  is  worthy  of  . llustraticn , see  Figures  36  and  3S.  It  is 
strongly  suspected  that  resonance  phenomena  involving  coupling  between 
turoulent  vortices  and  the  superposed  oscillatory  signal  is  present. 

In  [27],  and  as  seen  in  Figures  36  and  36,  the  experimental  data 
corroborates  the  use  of  a quasi-steady  model  for  dimensionless  frequency 
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rotics  iel2w  and  ur  to  the  transiticr.  region,  and  tr.e  use  of  a lar.ir.ar 
oscillator'.-  Docel  at  ratios  aoove  this  region.  The  exterirental  e-zi- 


cenoe  i.t  [2T,3i]  further  e.-'phasizes  the  futility 


develot  a rigorous,  ocntlicated,  sr.cotr.  transition  n;odel  as  -uas  cone 
in  attenuation  and  phase  -.relooify  behavior  in  the  transition 


rerocr.  is  tar 


fror.  smooth.  Fortunately,  the  transition  frequency 


bandwidth  is  predictarle  [27],  and  can  be  avoided  in  any  attempted 
transmission  cf  information  t.nrcugh  the  developed  turbulent  m;ean  flow. 
In  view  of  the  experimental  evidence,  it  is  thought  s'ufficient  to  use 
tne  constant  LRC  miooel  for  dimensionless  frequency  ratios  in  the  mid- 
range and  ut  tc  those  at  which  the  laminar  oscillator^’  model  becomes 


High  Frequency  Domain  - Turbulent  Mean  Flow.  In  the  high  frequency 
dom.ain  (w  >>  the  flew  becomies  frequency  dependent.  At  high 

frequency,  the  laminar  oscillatory  bcundar-y  layer  becomes  considerably 
thinner  than  the  equivalent  mean  flow  boundary  layer  <<  6v^)  3s 
previously  discussed  in  this  section.  Tnus  the  high  frequency  behavior 
becomes  essentially  independent  of  the  mean  flow,  and  approaches  the 
laminar  oscillator^’  behavior  associated  with  the  line  being  investigated. 
The  normalized  lam.inar  oscillatory  behavior  cf  both  circular  and 
noncirc’ular  lines  v.-as  shown  to  be  nearly  identical  when  plotted  against 
the  laminar  characteristic  frequency  ratio  u>v£  = --  reTr.ains  to 

determdne  the  transition  from  middle  to  high  frequency  and  to  normalize 
the  high  frequency  domain  results  such  that  they  are  applicable  to 
circular  and  noncircular  lines  carrm'ing  laminar  or  turbulent  mean  flow 
over  a practical  range  of  Reynolds  number.  Results  are  sought  which 
are  similar  tc  the  single  cur-<^es  of  Figure  37.  These  curves  include 


.15 


bcth  lov  ar.i  ::.id-freq uer.cy  perf orr-ance  of  circ-iar  ar.d  ricrjcircular 
liTi'sE  cauylr.g  = ccTr-prassih-e  or  ir.co-.pressidle  flui:  a~  t:,e  specified 
reynclds  .'UTd.er. 

If  Kichols'  ecuctior.s  (l‘»2  through  1*^1)  are  used  tegsther  with  rv » 
the  predicted  high  fre:y-er.cy  attenuaticri  is  toe  high,  as  showr.  ir. 

Figure  jI.  Thus  r^_  car.nct  be  used  as  the  characteristic  radius  in  the 
high  frequency  region  unless  ether  adjustaents  are  ir.ade.  Experimental 
evicence  supports  the  use  of  the  lar.ir.ar  oscillatory  model  in  the  high 
frequency  regions,  and  the  following  development  will  provide  the 
mathematical  relationships  and  adjustments  needed  in  the  high  frequency 
domain. 

Equation  (!“*<)  may  be  expanded  using  equations  (1^*2)  and  (1^3) 
with  V = vj  to  give 

1/2 


wnere 


(166) 


anfl 


(169) 


The  normalized  impedance  ratio  is  obtained  from  equations  (142),  (143), 
(145),  and  (36); 


zo  r 1 r 1 

= 1 - J(^)  1 + (>"-1)  J(4,)  = h(4i) 

^so  L J L J 


(170) 
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variation  of  Unad  justed  Fretiuency  Pepr>n(l«»nt  Attenuation  witli  Roynold.':  Nuniliof  in  t'irrnl.n 
and  Noncirr.ular  Lines 


',r- 

1 


w 

i 

1 


It  is  clear  that  the  arg'ccier.t  t cor.pletely  deterr.ir.es  t.te  hracheted 
extressior.  cf  equation  (161)  and  all  of  equaticn  (1?1).  The  r.odulus 
of  « car.  re  revritter.  in  terris  cf  the  characteri s tic  radius  and 
frequency  as 


(171) 


where  r.y  = rc,,  r- 

or  ry.^,  and  sim. 

ilarly  for 

Ljy . The  argument 

$ depends 

only  on  u/uy,  or 

its  equivalents. 

However , 

from  equation  (167 

) it  is 

seen  that  f also 

depends  on  a:/c^. 

Supposing 

two  lines  cf  the 

same 

characteristic  frequency  ratic,  but  different  wy,  i.e., 


i£_ 

U)\; 


1 


and 


WVj  r 


Lies 


0)2  wv2 


, for  ail  u. 


(172) 


Consequently,  a (and  £)  may  be  adjusted  between  lines  of  different 
characteristic  frequency  tc  obtain  the  proper  values; 


It  follows  that  a single  baseline  curve  can  be  constructed  to  describe 
the  high  frequency  attenuation  of  all  lines,  carrying  either  laminar 
or  turbulent  mean  flow.  One  fcnr.  of  this  curve  is  obtained  by 
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norrr.aliiir.g  a,  equation  (167),  in  a nianr.er  sir.ilar  to  Brown  [1:]: 


when  a = a(r^,).  Equation  (lEu)  shows  that  the  ncrr.alized  attenuation  is 
a function  of  only.  Equation  ''I?-)  is  plotted  in  Eig'jre  i-l.  The 

high  frequency  behavior  as  shown  in  Eig’ure  41  (above  u)/u,„,  = 2),  nay  be 
adjusted  to  the  appropriate  lanir.ar  reference  (aswiptote)  by  using  the 
knovT.  high  frequency  behavior  arid  the  approach  leading  tc  equations  (1"2) 
and  (173),  or  by  noting  that  the  correct  characteristic  frequency  ratio 
at  high  frequencies  is  whether  the  underlying  steady  flow  is 

larrinar  or  turbulent,  and  interpreting  figure  accordingly. 

It  is  convenient  to  norrr.alice  a such  that  the  onset  of  frequency 
dependent  behavior  begins  at  unity  value  on  the  Tr.id-fr'ecuency  plateau 
and  joins  the  laninar  oscillatory  liir.it  curve  of  the  fluid  line  being 
considered.  The  nomialited  high  frequency  attenuation  is  to  be  detenr.ine 
in  terms  of  The  appropriate  asymiptotes  are  given  in  equations 

(157)  and  (96)  for  the  middle  and  high  frequency  regions,  respectively. 
These  are  used  tc  define  a new,  extremely  simple,  single  expression  for 
the  middle  and  high  frequency  dimensionless  attenuation.  The  two 
asymptotes  are: 


“‘K/o'K,  = * 


(175) 


(175) 
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Tr.is  furiC'ior.  r<5preser.ts  the  rr.iddle  ar.d  high  f 
atter.-^ation  of  lar.inar  oscillatory  signals  in 
lines,  carrying  fully  cevelcrec  lar.inar  cr  tur 
inclusion  of  y in  this  equation  also  accounts 
Equation  (176)  is  plotted  in  Figure  u2.  The  h 


recuencv  direns denies s 
circular  and  ncncircular 
•huler.t  thrcughflow.  The 
for  coiT.pressihility  effects, 
ligh  corner  (break) 


frequency  occurs  at  the  intersection  cf  the  two  asymptotes--. as  previously 
shown  in  the  development  cf  equations  (161)  t.nrough  (165).  A new 
dimensionless  grouping,  the  abscissa  of  Figure  -2,  is  formed  which 
combines  mean  flow,  geometr\’,  and  frequency  dependent  paranet«rs.  This 
new  grouping  facilitates  comparisons  with  previous  work,  notably  that  cf 
Funk  and  Wood  [If],  as  discussed  in  the  following  s’ubsection . 


Comparison  cf  High  Frequency  Results  with  Previous  Work 

In  comparing  the  above  results  with  those  of  Brown,  £t  a^.  [15,27], 
and  Funk  and  Wood  [25,26],  it  is  necessary  to  establish  a common  basis 
among  the  references.  This  basis  assum.es:  (1)  lines  of  circular 
cross  section  (fvj  * r^  or  R);  (2)  nearly  incompressible  fluids  (y  = 1); 
(3)  unity  Prandtl  Number  (o  = 1);  (4)  constant  phase  velocity  (c  = c^); 
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(:)  develoz.ed  turbulent  -lear,  flow;  ( c)  relatively  snicll  airclitude 
ostillatcr"-  signals;  C’)  snail  attenuaticn  ter  line  wavelength;  and 
(£)  hign  ;r  near  high  signal  frequencies  > 1),  for  all  turbulent 

r.evnclds  numbers.  Since  runh  and  Wood  [2  5]  have  cent  area  their  work  tc 
tnat  of  crown,  £2  il'  l15],  the  comparison  herein  deals  primarily  with 
reference  [2:].  The  dimensionless  paraneters  of  Brown,  equations  (152) 
and  (153),  will  be  used  for  overall  comparisenE , since  these  parairieters 
nave  found  wide  use  in  the  literature. 

Tr.e  model  developed  by  funk  and  Wood  [25]  assigns  all  viscous 
effects  to  a single  fictitious  lan.inar  bcundarn’  layer  whose  thieVnaess 
(1)  is  based  or.  steady  flow  laminar  and  turbulent  friction  factor  data: 


d 16 


(1T5) 


where  R is  the  radius  of  the  circular  line.  All  other  effects  are  assigned 
to  an  inertial  plug-flow  cere  (i.e.,  L - La).  The  model  makes  use  of 
the  eight  assumptions  listed  above,  and  results  in  a single  dissipation 
function  (?)  such  that  the  attenuation  can  be  expressed  as 


a 


V? 

RCgb 


(180) 


The  attenuation  itself  is  assumed  to  be  small,  and  only  first  order  effect 
are  considered  in  the  expansion  of  the  dynam.ic  pressure  ratio: 


b?. 


-02  , 
e =1-02 


(181) 
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T>6  rrocel  of  rur>.  and  Wood  still  reouires  the  cotipotctior.  of  a cotitlex 
trar.sfer  fuoctior.  relatir.e  the  local  volocitv  gradier.t  to  the  cere 
velocitv.  ar.d  also  requires  r.'ut.erioal  ir.tegraticr.  of  the  dissipatior. 
fuTiCtior.  across  the  boundary  layer.  3o.h  of  these  requirer.er.ts  are 
avoided  by  the  tiodel  given  in  equation  (17S)  and  plotted  in  figure  u2. 
Fret,  equations  (65),  (75),  (7?),  and  (lUl)  it  is  easily  shown  that  for 
the  sate  circular  line  carrp'ing  turbulent  mean  flow. 


“vj.  _ _ (fPe)^  _ 

^ ^ r.  ' ^ r-v^  ^ (f?.e)t  ""  (rT.ey).t 


(1£2') 


The  nodel  of  Funk  and  Wood  and  the  author's  model  are  related  tnrough 
the  turbulent  mean  flow  parameters.  Comparing  equations  (179)  and  (1S2), 


i 1 / * vt  s' 

R “ r V ~F~  / 


(183) 


The  relationship  between  i and  the  para.meters  developed  herein  is  obtained 
by  comparing  expressions  for  the  dimensionless  attenuation.  From, 
equation  (180), 


^ ^a  R . 

— a = ■■  ■ • = — ft 

V V RCaii  A 


(16-4) 


and  from  equations  (176)  and  (163), 


awau^  = 4 


i^)  ““t  = I “H, 


(165) 


It  follows  that  equivalence  is  of  the  form,  i - 5^^.  From,  equations 
(178),  (183),  and  (185),  with  y = a = 1,  and  ryj^  = F-  for  comparison 
with  [25], 
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Lir. 


Lir.  r 

tL-«  ii 


(Ic?) 


and  at  Low  frecuency  ratios, 


Lin 

ur*-0 


R 

A 


(186) 


arid  The  high  oorrier  (break)  frequency  occurs  at  L - Fron.  Nichols  [3] 

and  Karar.  [20],  the  high  frequency  lariinar  oscillatory  attenuation  was 
approximated  by  equation  (96).  Equation  (96)  reduces  to  equation  (187) 
as  shown  in  the  following  development  (y  = o = 1). 


t 


y-1 

o 


m 


/ 


(189) 


Using  the  definitions  of  a>v£  and'6^^,  together  with  equation  (169)  and 
Brown's  notation,  gives 


r2c. 


V 


“H 


(190) 
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It  rerr.ains  tc  shew  the  ecuivaler.ce 


f the  ~id-fi'etuer.cv  asvmttctes , 
c = 1: 


The  last  twe  terr.s  of  the  areve  equaticr.  aa^e  ider.tital  tc  those  of 
eo’aatior.  (1£3). 

The  attenuation  Tr.odel  err.hodiec  in  equations  (ITt)  end  (165), 


atplicatie  tc  the  Tr.iddle  and  high  frequency  regions,  constitutes  a verv 
significant  sintlificaticr.  when  cct.pared  tc  the  models  given  in  [15]  and 


[25].  An  attreciaticn  of  this  model  can  only  be  realized  bv  comcarison 


with  the  cited  references. 

The  results  of  rrown,  et_  [If, 27],  Funk  and  Wood  [25],  and  the 
author  are  plotted  for  comparison  in  Figure  43  at  two  Reynolds  nuiribers. 
Re  = 1x10**  and  Ixio^,  for  r^^  = R.  In  view  of  the  experimental  data, 
the  author's  model  is  as  valid  as  any  of  the  referenced  models  for 
circular  lines,  and  is  comparatively  much  simpler.  Furthermore,  the 
model  is  applicable  to  any  fluid  line  of  arbitr'ary  cross  section, 
carrying  laminar  or  turbulent  mean  flow.  The  other  models  are  not. 


Turl)ijlenl  Mean  I'low 


VI. 


i.-.€  .ransier.t 


I I**-.' 


» 
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ii  r^^7-r£ 


i.'ev5j.cze-  iUrau.er.t  .-.ear.  : _ow 


In  this  section,  the  high,  rr.iidle,  end  Ic-k  fretuency  models  develtpet 
in  Section  V are  transferred  into  the  time  domain  using  the  general 
attroach  of  .Karam  [1P,2C].  The  methodology  and  form  of  the  data  is  verm- 
sim.iiar  to  that  presented  in  Section  IV;  much  of  the  general  discussion 
presented  in  Section  IV  concerning  the  various  regions  cf  the  frec.uency 
and  time  domains  is  alsc  applicahle  here,  /.cccrdingly , only  major 
differences  viil  be  pointed  out  in  the  development  and  discussion  which 
follows.  In  all  developtients  of  this  section,  unity  Frandtl  number  is 
assumed . 

Frequency  and  Time  Domain  Equations 

The  frequency  doir.ain  is  divided  into  three  regions  (low,  middle,  and 
high)  based  on  the  turbulent  characteristic  frequency  ratio,  - 

Ry^/La.  The  regions  are  divided  according  to  the  value  of  the  low  and 
high  corner  (break)  frequency  ratios,  as  defined  by  equations  (15C), 

(157),  and  (162): 

u: 

~ ' '“vt  I 


Low  frequency: 


uiv-*-  2 u). 


1 u;  4 “'^t 

Middle  frequency;  x < < 

•4  U) 


.4ign  frequency; 


U) ^ 4 ^ 


(192) 


y u, 


vji  "Vt 
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It  is  seer,  that  the  rericr.s  are  r.izhly 
arc  t’-rrbuler.t  steacy  flow  tarao.eters. 

vised  tc  recefir.e  the  direr.si 


deter.detit  or.  hcth  the  la-tinar 
A r.ew  characteristic  tiir.e,  t^  = 
tr.less  tir.e  tarar.eters  , such  that 


W'  1 ^ 


Wvt' 


:/t. 


(193) 
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Th€  ahcv€  frecu^ncv  T'®r*lcns  snc  diTr.6risicr.l6ss  tirie  parsTT}6t6rs  cr6 
used  ir.  develctirig  expressicr.s  for  the  various  pressure  and  flow 
■tr&iis  16^^  r‘6SDcns6s  *‘r‘6s6n't6d  ir.  *this  S6c*ticri« 


Line  Attenuation  Factcr 

To  facilitate  the  use  cf  .Karair. 's  r.ethodology , the  attenuation  in 
turbulent  mean  flow  is  expressed  in  tenris  of  the  three  logarithir.ic  line 
segments  (asymptotes)  corresponding  to  the  three  frequency  regions  cf 
equations  (192): 


X If 
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“I 


- K * V (ij 


2 


(194) 


“H* 


Ky: 


It  is  seen  that  the  high  and  low  frequency  attenuation  are  of  the  form, 
a = K /w.  Consequently,  the  solutions  of  Metzger  and  Vabre  [23]  and 
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:.arar.  [:c3  r.ay  be  acarted  :c  obtair,  the  vari; 


•us 


ar.ciAr*  r 


resz-cnses . 


rhe 


Etzer.uaticri 


ci-'ciu*  is  redefined  3S 
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au 


rcr 


(195) 


r*ere 


Pressure  (Floy)  Respor.se  tc  a Step  Input  of 
Pressure  (Floy).  Turbulent  Kean  Flow 

T.te  general  solution  for  the  transient  downstreati  pressure  (flow) 
response  to  a step  input  of  pressure  (flow)  was  given  by  equation  '102). 
Using  the  attenuation  ratios  of  equations  (195)  above,  together  with 
equations  (102)  and  (1S3),  the  generalized  transient  response  may  be 
expressed  as 


= erfc 


V^vj, 


to 


/2(t-to ) 


U(t-To) 


(196) 


L -I 

where  subscript  t denotes  applicability  to  turbulent  mean  flow  conditions. 
The  dimensionless  step  responses  for  long  and  short  times  corresponding 
to  equations  (1C7)  and  (106)  are,  respectively: 
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= errc 


1 (bh 
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-*  \imcJ  '■ 


'J  ( 1 -T  r ) 


= errc 


')  ' 


T-Tq  / 


U (t-Tq) 


(127) 


and 


( 13£ ) 


Ter  cor.r  letenesE , the  intermediate  or  middle  time  solution  is 


erf  c 
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. m 
KT  To 


►^2(t-To) 


U (t-to) 


(199) 


It  is  seen  that  the  middle  and  long  time  transient  response  is  completely 
determined  for  a given  fluid,  T,and  tq.  It  is  also  noted  that  equation 
(197)  is  identical  to  equation  (1C7);  in  the  latter  equation  the  times 
are  nondimensionalized  with  respect  to  Wq.  Equations  (197)  and  (199) 
describe  the  transient  response  of  arbitrary  lines  filled  with  a 
compressible  fluid  carrying  leurinar  or  turbulent  mean  flow.  For  laminar 
mean  flow,  is  used  to  define  the  dimensionless  time  (i.e., 

T = u-vj^t);  for  turbulent  mean  flow,  = Rv^^/la  is  used.  The  short 
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•cir.t  respcr.se,  ecuetrcr.  ^i9£),  retair.s  scr.e  iepe 

the  turbulent  to  la.t.lr!ar  characteristic  racii.  

aticn  is  that  the  high  corner  creak  frequency,  which  essentially  r. 
,he  division  of  influence  he tween  steady  turbulent  near  flow  hehav' 
and  frequency  dependent  larr.inar  oscillatcry  behavior  must  include 
•■nfluences  from  both  phenomena.  The  ratio  essentially  ad- 

he  transient  response  to  account  for  the  various  arbitrain-  line 
gecmetries  and  t'urbulent  Reynolds  numbers,  and  is  therefore  net  sc 
dependent  on  the  turbulent  characteristic  frequency  (u-v-). 

Equation  (196)  is  plotted  in  Figure  as  a function  cf  i = u 
for  tv;c  values  of  the  ratio  which  correspond  to  the  folio 

abulated  conditions  for  circular  lines  (ry^^  = R). 


dependent  or.  the  turbulent  characte 


abulated  conditions  for  circular  lines  (ry  = R). 


TAELE  VI I 

Equivalent  Circular  Line  Steady  Flow  Parameters  for 
FiEures  u4.  45.  and  46 


Tyj 


.44975 


1.4017 


Re 

fRe 

J(— ) 

“ '''vt/ 

1x10** 

316.4 

.15760 

1x10® 

1779.2 

.02603 

1 

F.(=ry^) 


,07996  1.4017 


The  info-uence  of  on  the  short  time  response  is  smal^.  and  nearl 

negligible,  as  can  be  seen  in  Figure  4u.  in  the  middle  or  long  time 
regions,  there  is  no  influence,  which  is  evident  from  equations  (19?) 
and  (199).  For  the  large  differences  in  mean  flow  parameters  given  in 
Table  VII,  the  nearly  unnoticeable  effect  on  the  pressure  and  flow 
param'eters  plotted  in  Figure  44  indicates  that  the  transient  responses 
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Figure  . rreKr:ure  (Flow)  Response  to  d Step  Input  of  Fressiire  (ft. 

of  Arbitrary  Cross  Section,  Carrying  Turbuleni  Moan  Flow 


=re  V'ii'v  weak  furiCtions  cf  The  equa'lrns  iavel;,pea  ir.  tr.is 

secticr.  alac  apply  to  upstrear  anl  dcwTiStrear.  flow-  transients;  i.e.. 


'J(t)  = F£(t)u(t)  cr  Q5(t)u(t) 

as  appropriate.  Expressions  for  the  re.t, airing  -jnit  step  inputs  and 
responses  of  interest  are  developed  in  the  next  suhse-STrSn. 

Pressure  (Flow)  Response  to  a Step  Input  cf  Flow 
I Pressure ) , Turbulent  Kean  Flow 

The  two  rericining  ditiensionless  transfer  functions  cf  interest  are 


^sc  \ 5s(s  ) y s ZgQ 
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1 ( 1 -dPtCs) 

[f^rn  ■ S Yso  ^ 

t 


(2C1) 


As  was  done  in  Section  IV,  the  approximations  used  by  Karam  [20 J are 
adapted,  witn  t = u'v.^t,  and  equations  (113)  and  (114)  are  used  to  obtain 
the  final  time  domain  expressions.  Subscript  t again  denotes  applica- 
bility to  turbulent  mean  flow  conditions. 

Approximations  of  the  characteristic  impedance  and  admittance  for 
the  long  time  (small  s)  and  short  time  (large  s)  regions,  adapted  to 
turbulent  mean  flow  conditions,  are  respectively: 
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As  before,  it  is  ccnvenierit  to  coinbirje  the  impedance  and  admittance  into 
a single  function  applicable  to  the  entire  time  interval: 
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( 206 ) 


Y^(X)dX  = 


(209) 


Usir;g  the  fcrrr.  cf  equaticns  (113)  and  (ll^^)  together  with  equations 
(196),  (2C6),  and  (209),  the  desired  dimensionless  step  responses  are 
obtained  as  fellows; 


(210 


(211) 


Equaticns  (210)  and  (211)  are  plotted  in  Figures  and  46,  respectively, 
for  the  paraJTieters  given  in  Table  VII. 

hysical  I.m;  li  cat  ions  cf  the  Fluid  line  Step 
esponse  with  Turbulent  Kean  Flow 

It  is  of  interest  to  discuss  the  physical  implications  of  the 
dimensionless  step  responses  shown  in  Figures  44,  45,  and  46.  For  lines 
with  the  same  value  cf  t^,  the  computed  responses  will  be  identical 
for  values  cf  t 2 (y/**)  (wv^/uiv.^ ) , and  will  be  nearly  identical  for 
smaller  values  of  t.  However,  two  lines  of  identical  tq  and  working 
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fluid,  but  differer,*  turbuler.t  character istic  frequencies  necessarily 
iraply  different  physical  lenjths  down  the  lines.  The  line  having  the 
larger  turluler.t  character istic  frequency  will  exhibit  a faster  rise 
to  a giver,  ratio  cf  ?r/^s  C'r/Qs?  a shorter  length,  and  will 
apt  roach  a diffusive  RC  line  ir.ore  rapidly.  Corritaring  the  two  circular 


^neuc.atic  lines  for  the  conaitions  of  Table  VII  with  t 


C - ‘ ‘ 


seen 


that  an  increase  in  Reynolds  number  from  Re  = 10“*  to  Re  = IC^  decreases 
tne  adiabatic  delay  time  (T  = i'c^)  by  the  ratio  of  31.6/1;  The 
implication  cf  this  example  is  that  increases  in  the  turbulent  flow 
Reynolds  nu-mier  greatly  accelerate  the  transition  to  diffusive  RC  line 
behavior;  consequently,  high  turbulent  Reynolds  numbers  should  be 
avoided  in  any  attem.pt  to  transmit  infcnnation. 

In  terms  of  the  equivalent  lam.inar  and  furbulent  steady  flow 
boundary  layers,  equation  (163),  it  is  seen  that  the  condition 


T 2 


1 


is  ecuj.vaj.ent  to 


T 2 y ^ 


for  circular  lines  or 


T 2 


^VC 


for  arbitrary  lines.  Thus  the  break  time  between  short  time  and  inter- 
m.ediate  tjm:e  domiain  responses  is  a strong  function  of  the  equivalent 
steady  flow  boundary  layer  parameters. 
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VI I . CcriClusi jr.s  and  Recorr.er.dations 

This  section  contains  a sun.marn’  of  the  conclusions  reached  as  a 
consequence  of  this  work,  r-ecoriiaendations  for  extensions  of  this  work 
are  addressed. 


Conclusions 

1.  In  general,  the  ver>’  close  analogic  between  the  STnall-signal 
response  of  fluid  arid  electrical  transirission  lines  has  been 
reinforced  for  a wide  range  of  fluid  line  gecr.etries , Ttean  flow 
conditions,  fluid  media,  and  heat  transfer  effects.  The  signifi- 
cant role  of  the  equivalent  electrical  trar.sr.ission  line  properties 
in  describing  the  fluid  line  response  has  been  demonstrated. 
Distributed  line  parameters  such  as  resistance,  inductance, 
conductance  and  capacitance,  derived  from,  the  series  impedance 

and  shunt  admitta.nce  of  an  equivalent  electrical  transmission  line, 
are  shown  to  be  the  key  paramieters  in  detenr.ining  the  line  response. 

2.  A closed  form.,  distributed  parameter,  analytical  solution  for  the 
small-signal  frequency  response  of  fluid  transmission  lines  of 
annular  cross  section  has  been  obtained.  The  solution  considers 
shear  flow  effects,  fluid  compressibility,  and  heat  transfer 
effects,  and  is  applicable  to  blocked  lines  or  lines  carrying 
developed  laminar  mean  flow.  The  solution  was  shown  to  be  in  very 
good  agreement  with  experimental  data  for  blocked  annular  pne-jmatic 
transmission  lines  over  a wide  range  of  annular  radius  ratios  and 
signal  frequencies.  In  comparing  the  annular  line  results  with 
those  previously  obtained  for  rectangular  lines,  it  was  shown  that 
for  large  radius  and  aspect  ratios  the  solutions  converge.  The 
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ar^nular  line  scluticr.  contains  tne  circular  line  sclutio 
special  case. 

3.  lines  of  arlitrary  cress  section  can  re  accuratelv  ttcceled  as 
equivalent  circular  lines  having  the  same  ratio  of  steady  flow 
resistance  to  inertia  properties.  The  concept  cf  an  equivalent 
characteristic  frequency  and  radius  can  be  applied  to  transmission 
lines  of  arbitrary  cross  section  to  predict  beth  the  frequency 
response  and  transient  behavior  of  these  lines,  including  either 
laminar  cr  turbulent  m.ean  flow  effects.  In  solving  practical 
problems  involving  fluid  system.s,  this  equivalence  should  be  of 
significant  value  in  simiplifying  the  system  miodeling  task. 

4.  The  characteristic  radius  is  analogous  to  the  hydraulic  radius 
used  to  represent  circular  and  noncircular  lines  in  steady  flow 
analysis.  However,  the  characteristic  radius  mjcre  accurately 
accounts  for  the  actual  viscous  and  inertial  effects  and  was  showTj 
to  be  clearly  superior  to  the  hydraulic  radius  in  modeling  the 
frequency  response  of  noncircular  lines.  In  view  of  the  much  closer 
agreement  with  existing  data  and  analytical  solutions  for  both 
rectangular  and  annular  lines,  the  characteristic  radius  should  be 
used  instead  of  the  hydraulic  radius  to  determine  the  equivalent 
circular  line  for  all  noncircular  lines . 

5.  Very  significant  simplifications  have  been  achieved  regarding  the 
modeling  of  the  frequency  and  temporal  response  of  circular  lines 
carrying  turbulent  mean  flow.  Comparison  of  the  simplified  models 
with  existing  models  and  experimental  data  on  circular  lines  shows 
that  the  simplified  models  are  as  accurate  as  any  existing  models 
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Ln  predictir^g  the  frequericy  response,  rurthertr.cre , the  si-plified 
rr.odels  are  net  restricted  tc  circulsr  lines  cr  inct-Trpreesihle 

fluids . 

£.  Of  the  four  equivalent  circuit  cort'cnents.  the  steady  flew  viscous 
resistance  is  the  most  influential  parair.eter  in  deterr.ining  tne 
frequency  and  teTnporal  respense  of  fluid  transrussicn  lines.  For 
developed  laminar  flow,  the  steady  flow  resistance  is  determined  by 
the  line  geometry  and  fluid,  and  is  not  influenced  by  changes  in 
the  F.eynclds  number.  For  developed  turbulent  irsan  flow,  tne 
resistance  is  strongly  dependent  on  the  Reynolds  n-umuer.  In  all 
cases,  the  thickness  of  the  frequency  dependent  laminar  oscillatory 
boundary  layer  relative  to  the  equivalent  steady  flow  boundary 
layer  is  a key  parameter  in  determining  whether  the  frequency  or 
transient  response  exhibits  predominantly  viscous,  iscthermial 
behavior  cr  inertial,  adiabatic  behavior. 

Recommendations 

1.  The  small  signal  transient  response  of  both  circular  and  non- 
circular  lines  carrying  laminar  throughflow  should  be  deterir.ined 
experimentally  and  compared  with  the  analytical  results  presented 
in  this  work.  The  experLment  should  include  both  compressible 

and  incompressible  fluid  flow,  and  should  be  repeated  for  turbulent 
throughflow . 

2.  More  analytical  and  experimental  work  is  also  needed  to  extend  the 
applicability  of  the  concepts.  The  frequency  and  temporal  responses 
of  fluid  transmission  lines  to  large  amplitude  excitation  signals 
are  of  particular  interest.  The  assumption  of  essentially  constant 
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Ar-Lendix  A 


DC  Pesistajice  otjC  DC  Inductarjce  of  7 
Lir.es  Carr'.'ir.r  DevrLcred  Lar 


-urc  Tr ar.srr.rssscr. 
ir.ar  riov: 


In  g^Tieral,  the  DC  resistance  (?.•,•)  and  DC  inductance  (Ly)  are 
cetem.ir.ed  free,  the  steady  flow  parar.eters  cf  the  rigid  trarisr.ission 
line.  These  parameters  will  be  functions  of  the  line  geometry,  fluid 
properties  and  steady  flow  profiles.  The  DC  resistance  is  defined  by 


-SP/Pz  -cF/Pz 


(A-1) 


and  the  DC  inductance  is  defined  as 


' i (t) 


(A-2) 


where  u is  the  steady  flow  velocity  (either  laminar  or  turbulent);  u 
is  the  average  velocity,  u = 1/A  judA;  and  = p/A  is  the  adiabatic 
inductance. 

For  circular  lines,  Ry  and  Ly  are  constants.  Per  steady  developed 
laminar  flow,  Ry  aiid  Ly  are  obtained  from  the  classical  Hagen-Fciseuille 
flow : 


Ky  - r 


-3P/&2  StTm 


A‘ 


(A-o) 


The  Kagen-Poiseuille  velocity  ratio  is  [3C] 
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(A-4) 
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where  r is  the  radial  cocrdir.ate  ar,d  is  the  inner  radius  of  the 
circular  line.  Using  this  ratio  in  equation  (A-i')» 


(A- 5) 


For  annular  and  rectangular  lines  and  Ly  will  depend  on  the  radius 
or  aspect  ratio,  respectively.  For  the  annular  line  [£], 


- /i 


(-dr/dz) 


(.<■.-£) 


Fror.  which  Rvj-^  is  easily  obtained  as 
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-er/oz 
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(A-7) 


where  r^'  = is  the  radius  ratio  of  the  annular  cross  section. 

Tne  velocity  ratio  is  [ 8 ]: 


u 


u 


[j__r»>2]  In  pi  _ [l-rV^]  In  r*- 
4 [(itrf‘)  In  rV  + (i-r-^)] 


(A-e) 


wnere  r*  = r/rQ  and  the  dependence  on  geczietry  is  obvious.  The  notation 
r*  = ^i/^o  i“  generally  used  in  the  body  of  this  work.  However,  for 
this  appendix  only,  it  is  necessary  to  distinguish  between  the 
cumensionless  ratio  r>'  = t/Tq  and  r;  = Using  equations  (A-2) 

and  (A- 6), 
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The  DC  resistaTiCe  for  the  rectangular  line  is  easily  found  as 


Tr.e  velocity  ratio  is  giver,  by 


:£/>;r/7— 1 


. r.6  ir.cuc*tance  i.s 
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(A- 15) 


‘.6  ir.tegr-arid  is  of  the  foroi 


I y£(n)x/(C) 


= .1.  I yi(h)xi(C)yj(n)xj(0  (A-16) 


where 


Xi(;)  = 


i 

cosh  — C 


a,- 


cosn 


ana 


sinaj 


cos  a 
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iOe  integrand  is  sufficiently  well-behaved  to  permit  interchanging  the 
orders  of  sataation  and  integration.  Thus,  the  integral  cf  equation 
(A-15)  may  be  written  as 


OD  OC 
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^ = I I \ yi(n)y4(n)  dn 
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Xi(C)Xj(C)d5  (A-17) 
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intefei atlor,  is  greatly  sir.plified  zy  r.otir.g  that 


y-(')y--(Od-  = i *■  : 


(A-IS) 


rhe  integral  is  evaluated,  ar,d  Iv,.  becomes: 


4d‘  i = I 0 ^ i 


* . - ^i  \ 

tann^  — j 


(A- 19) 


It  is  of  interest  to  include  the  case  cf  steady  developed  latinar  flow 
between  parallel  plates,  becaute  this  case  represents  Uniting  values 
of  both  radius  ratio  for  annular-  lines  and  aspect  ratio  for  rectangular 


lines . As 


r-'  -*■  1 (a  - 0),  the  pairs  Ry^ » Ly  and  Ry  , Lv_  will 

li  n ^ 


approach  the  values  obtained  for  parallel  plate  flow.  Taking  a unit 
width  cf  flow  between  parallel  plates  where  2b  is  the  distance  separating 
the  plates  [31]  and  y is  measured  fret  the  center  plane  between  the 
plates. 
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And  the  DC  resistarice  per  unit  width  is 
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range  of  aspect  cr  radijs  ratios  for  rectangular  or  annular  lines. 

With  inside  radius  r^  = C for  annular  lines,  K.i,_  = K- __  = ^4/3.  As 
r^  -►  r^,  -»  1.2.  Fcr  rectangular  lines,  has  its  ir.a>:in:ui..  at 

a = 1 where  K*  = 1.3735,  a.nd  aptroeches  Yr  = 1.2  as  a ■*  0. 

A sir.ilar  approach  rjay  be  used  to  obtain  the  low  frecuencj'  (DC) 
viscous  laminar  resiEtar,ces  and  inductances  for  other  noncircular  lines, 
provided  the  steady  flow  velocity  profile  is  )<nown. 


Apperidix  E 

High  Frecuer.cy  Aprrcxir/atioris  for  Ar.r.ular  Lirjes 

First  order  at'proxir.ations  of  the  systet:  variables  desoribir.g  the 
terfcrr.ar.oe  of  annular  fluid  transTr.isf ion  lines  at  high  frequency 
ratios  >>  1)  are  obtained  frcrr.  equivalent  circuit  considerations 

and  analysis  of  the  as}Tiirtctic  behavior  of  the  Eessel  function  solutions. 
In  general  at  high  frequency  ratios,  the  flow  is  approaching  inertial, 
adiabatic  behavior.  As  was  p'ointed  out  by  .hichols  [;],  the  velocity 
profile  is  characterized  by  a large  central  plug  wherein  inertial 
effects  are  predominant,  and  a viscous  skin  regicn  wherein  viscous  shear 
effects  are  predominant.  The  viscous  region  is  characterized  by  a 
skin  depth  which  is  proportional  to  / v'/u . As  frequency  (iii)  increases, 
i-y)  decreases.  Similarly  heat  transfer  effects  may  be  separated  into  a 
relatively  large  adiabatic  core  region  bounded  by  S nonadiabatic  region, 
the  latter  being  characterized  by  a thermal  boundart'  layer  thickness 
6j  which  is  proportional  to  >'v.p/u).  For  the  annular  line,  there  are  two 
bcuncer-y  layers  arid  their  influence  on  the  velocity  and  temperature 
profiles  var;,'  with  radius  ratio.  The  behavior  is  described  in  the 
velocity  and  temperature  profiles,  equations  (11)  arid  (12),  respectively. 

/ I 

In  general,  the  Bessel  function  JoUj^ 
the  sclution  at  the  outer  wall  (r  = Tq 

influential  at  the  inner  wall.  At  high  frequencies  the  following 
inequalities  apply: 
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where  fcr  JqjJ;  the  argur.ert  (r)  is  presuined  to  be  rj"' and 

for  K'jK]  the  argufier.t  is  pres'jnied  to  be  ^v'lJTv  For  the  energy 

equations  v is  replaced  by  vt. 

With  the  foregoing  approximations  the  parameters  appearing  in 
equations  (11)  and  (12)  may  be  greatly  simplified.  The  coefficients 
£,  Z,  and  K of  equations  (13),  (14),  and  (16)  become: 
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From  [5,6]  the  first  order  approx imiat ions  of  the  Bessel  functions  are 
given  in  Table  B-I.  Using  the  expressions  given  in  Table  E-I  and  keeping 
only  first  order  terms,  the  high  frequency  series  impedance  becomies 


= j‘“La  jl  + 
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(E-6) 


and  the  shunt  admittance  is 


'ki  = l“Ca  jl  - (^)  ^ * 3(r-Hl)]j  (8-7) 
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I.eriv£ticr.  cf  Stt.cII  £i?r.=l  Ecuaticns 


sever; 


Cylindrical  Cocrcir.ates 

:r.is  appendix,  the  centinuity , eiorentur  and  energy  ecyaatiens 
: the  sr.d..  signal  prepagatien  cf  pressure  waves  threugh  a 
fluid-filled  line  are  developed  fror.  the  iTicre  general  ecyuations  cf  fluid 
mechanics  anc  theme cynami cs . The  ecyuaticr.s  are  expressed  in  a cylindric 
ccordinate  system  (r,  6,  z).  The  racial,  angular  and  axial  velocities 
in  tnis  cocrciinate  system  are  uj-,  ut,  ar.d  u-,  respectively,  and 
vector  velocity.  In  the  following  developments 
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The  principal  direction  cf  sir.all  signal  wave  motion  is  the  axial 
direction;  thus,  it  is  assumed  that  the  flow  is  axially  symmetrical.  Thi 
further  implies  Ug  = o/Se  = 0.  This  assumption  eliir.inates  the  momentum 
considerations  of  Equation  (C-2).  In  addition,  the  body  forces  (Fj,,  F^, 
F2,)  are  also  assumed  negligible.  The  momentum,  equations  are  thus  reduced 
to  the  following: 
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The  wavelengths  of  interest  are  assumed  large  compared  tc  the  cross- 
sectional  dimensions  of  the  line,  which  implies  that  the  radial  component 
cf  velocity,  u^,  and  its  derivatives  miay  be  neglected.  The  abcve 
e:, cations  are  further  reduced  to  the  following: 
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(C-£) 


This  assu-.iticr.  raav  be  expressed  as  the  following  condition:  /.  >> 
where  is  the  principal  dimension  of  the  cross  section  in  the  radial 
direction.  Since  c = fx,  and  f = a)/2r  this  condition  is  equivalently 
excressed  as 


G 


(c-e) 


where  c = speed  of  sound  in  the  line.  It  is  next  assumed  that  the 
pressure  is  uniform  across  any  line  cross  section.  This  assumtticn 
removes  the  dependence  of  F on  r and  eliminates  the  radial  momentum, 
equation 


(C-9) 


It  is  necessary  to  define  the  conditions  under  which  the  axial  variation 
of  Uz  is  negligible,  and  to  remove  this  dependence.  The  pressure, 
velocity,  and  density  m.ay  be  written  in  term.s  of  a mean  value  and  a 
smiall  perturbation  about  this  value  as  follows: 


C = p t p ' 

P = p + p'  (C-IC) 
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and  nhe  suiscritt  z has  been  dropped  frcr.  velocit\’,  since  only  the  axial 
velocity  remains  in  the  miOirientuii  equation.  Suhstituting  these  variables 
into  equation  (C-~),  and  separating  the  steady  and  dynair.ic  terns  gives: 
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Equation  (C-iu)  can  be  further  reduced  by  neglecting  higher  order  terms 
involving  products  of  perturbation  quantities; 
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(C-15) 


It  is  now  possible  to  solve  the  steady  flow  problem,  equation  (C-13), 
by  standard  techniques.  The  combined  homogeneous  and  particular  solution 
for  velocity  is  as  follows: 
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(C-l£) 


where  C;  and  C-  are  arbitrary,’  constants,  s'oiject  tc  the  boundary; 
conditions  of  either  circular  or  annular  lines.  For  circular  lines,  the 
boundary'  ccnditicns  are: 


u remains  finite  at  r = 0 
u = 0 at  r = r^ 


(C-17) 


Applying  these  conditions  to  equation  (C-16)  leads  to  the  well-known 
Hagen-Poiseville  velocity  distribution: 


u = - 
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3z  4u 


-rr  (rc^-r^) 


For  the  annular  line,  the  boundary  conditions  are: 

} 

u = 0 at  r = r.;  ( or  r*  = rp* ) 
u = C at  r = r©  (or  r*  = 1) 


(C-18) 


(C-19) 


where  r*  = t/Tq  and  rp*  = ^i’  ^o  inside  and  outside 

radii  which  bound  the  annular  cross  section  of  the  line.  Applying  these 
conditions  tc  equation  (C-16)  gives  the  steady  developed  laminar  flow 
velocity  profile  in  an  annular  region. 
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Ret\.rr.ir*g  to  the  cynair.ic  prohlez,  ecuaticr.  (C-15),  it  is  r.etessar-_;  to 
develop  the  cor.ditions  for  which 
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This  is  done  by  assuming  the  dynair.ic  pressure  to  re  of  the  following 
fcr~: 


, , -Tz  luit 

p'  = ipe  e- 


(C- 


2) 


arid  r is  the  propagation  operator  aefined  in  Section  II,  which  defines 
the  lengthwise  dependence  of  the  dynan.ic  pressure.  Lifferentiating 
equation  (C-22)  gives 
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To  a first  approximation , the  axial  sensitivity  of  the  perturbation 
velocity  is 
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(C-2U) 


Substituting  equation  (C-22)  into  (C-2u),  then  differentiating 


equation  (2-2^^)  with  respect  tc  z gives 
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(C-25) 
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;orr.pcring  equations  (C-21),  (C-23)  and  (0-25)  it  is  seen  that 
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(C-2c) 


For  air  at  standard  conditions  this  amounts  to  requiring  that 
u <<  3.2’'IC'  rad/sec.  Tne  radian  frequencies  considered  in  this  work 
are  well  witr.in  this  restriction.  Using  equation  (C-26)  in  equation  (C-15), 
and  rearranging  terms 
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Equation  (C-27)  is  identical  with  equation  (2),  Section  II.  The 
continuity  equation  is  considered  in  the  following  development. 

Continuity  Equation . Tne  continuity  equation  for  a general 
compressible  fluid  expressed  in  cylindrical  coordinates  is: 
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By  previous  assumption,  Uj-  = ug  = c/S6  = C;  thus 
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Expanding  equation  (C-29)  gives 
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Energy’  Ecuaticn . Ecr  a corcr-sssible  fluid  with  ccr.stant  therTTiai 
cor*Guctivit>s  k,  the  energy-  equation  is 


Dh  E?  _ 5q 
Dt  ” Dt  ~ Jt 


+ 4 + kV*'?  - V 


(C-31) 


•where  h is  the  specific  er.thalp}’,  q is  the  internal  heat  generation 
ether  than  by  viscous  dissipation,  i is  the  mechanical  or  viscous 
Gissipation  function,  and  is  the  radiation  heat  flux  vector.  For 
small  signal  propagation  the  dissipation  4-  is  ass'cmed  negligible. 
Radiation  effects  (q^)  and  internal  heat  generation  by  chemical  or 
other  like  mieans  (q)  are  not  considered.  With  these  assumptions, 
equation  (C-31)  reduces  to 
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Expariding  equation  (C-32)  in  cylindrical  coordinates. 
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iirice 


~T.e  axial  velocit\'  reir.air.s,  the  surscrirt  z is  cr-csysc  at 


this  pcirt.  As  before,  the  variables  are  defined  as  ir.ean  values  and 
sr,£ll  perturbation  cuantities: 
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These  variables  are  substituted  into  equation  (C-3'^}  to  obtain  the 
small  signal  energ\’  equation.  The  result  is 
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Using  the  relation  dh  = CpdT  and  dividing  by  the  quantity  dC-  gives 
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Tne  Frandti  nuirjb'er  = uC^/k,  together  with  the  relation  z/Z  = v,  gives 
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ar.d  vp  is  the  thermal  diffusivity.  Substituting  equation  (C-3c)  into 
(C-37),  and  rearranging 
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It  is  necessary  to  detem.ine  the  conditions  for  which 
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Assur.ing  the  dvr.ar.ic  terturbatior;  cf  temterature  is  cf  the  fol^ 
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Since  the  ir.oduius  of  r is  approximately  equal  to  id/  it  is  seer,  by 
comparing  equations  (C-42)  and  (C-43)  that  the  condition  cf 
equation  (C-UO) is  satisfied  when 
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For  air,  - .7  so  that  the  condition  expressed  in  equation  (C-i-u) 
is  nearly  identical  to  that  of  equation  (C-26).  The  final  form  of 
the  energN-  equation  is 
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whicr.  is  in  agreement  with  equation  (3),  Section  II. 
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.er.pora^  .-erivativc  cr  Der.sitv 


thc'ug.-.t  necessary  to  shov.'  the  derivation  cf  ec’caticn  (22), 
sirice  This  6cucTicr.  is  cr-cicl  To  Ths  flsvelcTir.sr.T  cf  ohs  lins  shur:T 
properties.  The  density  is  assunied  to  depend  on  pressure  and 
ten-,perature , ? = c(?,T). 


QC 


= f||^  aptfn:)  dT 

VC.-;,  ^ 


(C-^5) 


rhe  isotherrial  coircressibi lit\'  e is  defined  as 


E = 
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(C-,7) 


where  V ss  the  specific  voluirie.  The  coefficient  cf  thermal  expansion 


6^  is  defined  as 


= V (st)  = ■ p (stI 


(C--^c) 


For  an  ideal  gas,  £pT  = 1 and  eP  = 1. 

Ecuation  (C-E3)  may  be  rewritten  in  terms  cf  the  temporal 
derivative  as 


1=  ^ If  - = It] 
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From  the  definition  of  the  specific  enthalpy,  h: 


dh  = de  * d(E) 


(C-50) 


where  e is  the  specific  internal  energy.  For  a perfect  gas 
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dh  = C-dT  ar.d  de  = C,.d7 


(C-51) 


Using  equations  (C-51)  in  equation  (C-50)  gives 
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Rearranging  terir.s  and  using  equation  (c-ug)  together  with  the  ideal 
gas  relations,  £^7  = eF  = 1, 
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The  pressure  tern;  in  equation  (c-49)  niay  now  be  expanded  as  follows 
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Substituting  equation  (C-55)  into  equation  (C-49)  leads  directly  to 
equation  (22)  repeated  here  for  convenience: 
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The  development  of  the  remaining  equations  in  Section  II  is  relatively 
streightforvard. 
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